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: Elementarer Beweis für die Eindeutigkeit der Strömung 
4 in der laminaren Grenzschicht zur Potentialströmung 
E U=u,x2” mit m2() bei Absaugen und Ausblasen 

Von Rudolf Iglisch, Braunschweig 


‚Man beherrscht die Strömung in derlaminaren Grenzschicht um einen ebenen unendlichausgedehnten Keilvom 
7 Öffnungswinkel Br, zu der die Potentialströmung U = u, xm mit m = B](2 — ß) gehört, wenn man die ge- 
2 wöhnliche Differentialgleichung f’’’ + ff” + B(1— f’?) = Ounter den Randbedingungen f(0) = 0, (0) = 0, 
ie f(») =1löst. Daß dieses Problem für B O0 nur eine einzige Lösung besitzen kann, unter der Zusatz- 
bedingung f’(n) > 0 für n > 0 (das bedeutet eine Strömung ohne Rückfließen), wird in der folgenden Note auf 
elementarem Wege bewiesen. 


Given an infinitely extended wedge-shaped body with the flare angle ß - x, and the potential flow U = u, am 
round the wedge, the corresponding flow in the laminar boundary layer is known, if the ordinary differential 
equation F’’ + f’+B(1— f’?) = 0 with the boundary conditions f(0) = C, f’(0) = 0, f(w)= 1 is solved. It 
is proved, by elementary means, that the problem has a unique solution for ß > 0, provided that f’(n) > 0 for 
n > (that means flow without back stream). 


Le couant dans la couche laminaire de limite autour un coin plat, etendu & l’infini, de l’angle d’owverture Br, 
auquel correspond le courant potentiel U = u, xm avec m = B/(2 — P), est connu, si l’on resoud l’Equation 
differentielle ordinaire f’’ + ff” + B(1— f?) = 0 aux conditions marginales f(0) = C, (0) = 0, f(») =1. 
Que ce probleme pour ß > O ne peut avoir qu’ume seule solution ü la condition additionnelle f’(n) > 0 pour 
n > (celqui signifie un courant sans retour), est prouv& dans la note suivante dans une maniere Elementaire. 


TeyeHne B JaMUHAPHOM IOTPAHHYHOM CIIOE OKOJIO ILIOCKOTO ÖeCKOHEYHO HPOTSASKEHHOTO KIIHHA C 
yIIoM packpsTu4 P-n, KOTOPOMy COOTBETCTByeT NHOTEeHIMAAIBHOB TeyeHme U = u, zm 
ce m=ß/(2— ß), ocBauBaeTca NyTeM pemeHnn IPOCToroO AHudhepeHmmansHorO YPaBHeHHA 
ft” +14” + BL — f?)= 0 upu kpaeBuXx ycıoBuax f(0) = (C, f(0)=0, f(®)=1. B mpenuaraemoli 
3aMETKEe BlIeMEHTAPHEIM O0PA30M NOKA3EIBAETCH, UTO PACCMATPHuBaeMaA MPoÖNeMa B cAy4ae 
ß > 0 Mo:keT UMeTb AHIMb ONHO eNHHCTBEHHOE peImeHHe IIPH NONOJIHHTEJIBHOM yCJIOBHH fmM)>° 
aaa n>Ü (T.e., IOTOK Öe3 O0PATHO!O TeyeHnn). 

1. Einleitung 


In einer früheren Arbeit [1] !) wurde bewiesen, daß das Randwertproblem 

a ee. DE N 

ET Eee) 

für alle Parameterpaare (C, $) mit ö > 0 mindestens eine Lösung besitzt, die die Zusatzbedingung 
Ve ey < oo 2 ARE) 

erfüllt. In dieser Arbeit soll gezeigt werden, daß das Problem (1), (2), (3) für die angegebenen 


Parameterwerte höchstens eine Lösung f(n) besitzen kann; es ist also eindeutig lösbar. — Mit 
; dE SR 
Ed=f(n)s Er Bused alsE naachee url. ı 1.55 4) 
erweist sich das Problem gleichwertig dem folgenden: 
1.0 SV RER (5) Y 
F(6y=0; EEE WRITER a IN ET TREE 
DES FIN für IE I BSEFO FB FIIR), 


so daß nur für das letztere die Eindeutigkeit der Lösung sichergestellt zu werden braucht. Der 
Eindeutigkeitsbeweis wird auch den Parameterwert $ — 0 mit erfassen, der beim Existenzbeweis 
beiseite gelassen werden mußte (Vgl. [2]). 


2. Hilfssatz 1 
Für zwei Lösungen F,(f) «= 1, 2) von (5), die durch die Anfangsbedingungen 
Bier Et us.» 8. «7 l) 


1) Der Leser findet dort alles Notwendige zur Einführung in das vorliegende Problem. 
30 
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bestimmt sind, sind die Differenz der Umkehrfunktionen 


KA BAN Gr 1A 


sowie die mit # multiplizierte Differenz der Richtungstangenten bei gleichem F 
ANA NEE UP re (10) 
bei von & bis 1 wachsendem F monoton wachsende Funktionen von F, solange beide F, > 0 sind, 
— sofern mL 
Dei Inst, und) 2 rare (11) 


ist. (In mindestens einer der beiden letzten Ungleichungen muß das Ungleichheits-Zeichen 
stehen, da sonst die Lösungen identisch wären). 


Folgerung: F,(f) und F,(f) können sich für& <F < 1 nicht schneiden, solange beide 
Funktionen monoton wachsen; außerdem gilt dann für <F<1 


N EEE (12). 
Beweis: Für die Funktion 
U) EFineind eb ZI en (13) 
erhält man aus (5) 
dU 1— FF? 
Eee u) Tepe, gelburangeanrs ad ‚naht St: bei (14), 
daraus für 
Mi 
D(F) = U,(F) — U,(F) = Er (Ei Böen a ee (15) 
dD pa — F?) 
reelle Tre rrellinerulh) =) ka Erle (16). 
Stegen dD(«) { 
un ist infolge (8), (11) und (13) D(a«) >0 und IpT >0, falls 8 >0O oder a, — fi >0 ist; 
: 0 dD 
im Falle = 0 und f, — fi = 0 gilt jedenfalls sn 2 >0 für genügend kleines e. So- 


lange mit wachsendem F die Größe (16) positiv ist, wächst D(F) monoton mit F an, damit nach 


(15) auch F?— F}. Bei gleichem F ist dann stets die Neigung von F? gegen die f-Achse größer 
als die von F?2; also ist auch bei konstantem F die Differenz (9) eine monoton wachsende Funktion 


dD 
von F, so daß nach (16) IF >0für alle # zwischen a und 1 gelten muß. Somit ist der Aus- 


druck (9) für alle diese F-Werte eine monoton wachsende Funktion. Die analoge Aussage für (10) 
bedeutet nur das zwischendurch erhaltene monotone Anwachsen von (15). 


3. Hilfssatz 2 
Für jede Lösung des Problems (5), (6), (7) gilt für alle f > C 


FSU MIFSEN OF ET SP 


Beweis: Im anderen Falle müßte es einen /-Wert geben, an dem F = O ist. Dort vereinfacht 
sich (5) zu 

PP = DPI PFEFEN (18), 
falls ö > 0 ist; es könnte also höchstens ein echtes Maximum vorliegen, dem dann aber infolge 
der zweiten Bedingung (6) ein Minimum folgen müßte, was wegen (18) für den entsprechenden 
[-Wert unmöglich ist. — Im Falle ß — 0 kann nie F = 0 sein, da sonst F = const die zugehörige 
Lösung von (5) wäre. 

4. Hilissatz 3 
Solange f > 0 ist, gilt für jede Lösung von (5), (6), (7) 


BUD » aus 
damit erst recht RR een. Sl) 


RR P-<N 2 ge. PO er 
Beweis: Unter Berücksichtigung von (17) folgt (19) sofort aus (5). 


DA and 1 5 dh un ShA 5 a u Aal ala ll a 2 0 2 nn 
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5. Hilfssatz 4 
Es gibt keine zwei Lösungen F,(f) (v = 1, 2) von (5) unter der Nebenbedingung (7), für die 


Ri aast mit}, 07 RR AN a (21) 
und gleichzeitig 


Beweis: Mit der Abkürzung 
Wi,= EINE, nn (23) 
liefert die Subtraktion der beiden für die F, angeschriebenen Gleichungen (5) 


FR, — Ph, + P,M--FM+ RR,—rFF)—pR—P)=0, 
FAR, — PR) + RP — 9) + FAR,— FR) + Fi —F) 
+ IF, —F)) LIEB, - F) BR m) 0, 


FW HLFR, HF) + HF IW HARD) R+F)+AHi)W=0 N. 
Da infolge (21) und (22) W(f,) = 0 = W() ist, muß es mindestens einen positiven f-Wert geben, 
an dem W(f) z.B. ein positives Maximum hat. Dort wire W >0,W=0,W <0. Das liefert 
aber in (24) einen Vorzeichen-Widerspruch, da der Faktor von W negativ ist, wie man sieht, wenn 
man die geschweifte Klammer mittels (5) so umrechnet: 
ey: = —B—F(P,+F)W +, — PF}F, 
und (20) beachtet. 
6. Der Eindeutigkeits-Satz 
Das Problem (5), (6), (7) hat bei # > 0 höchstens eine Lösung. 
Beweis: ImFallO > 0 folgt dies sofort aus Hilfssatz 4. Für © < 0 — die folgende Betrach- 
tung gilt natürlich auch für CO > 0 — betrachte man eine Horizontale #=cmitl— e<ce<1 
bei genügend kleinem e. Unter Annahme zweier Lösungen F, und F, wird diese Horizontale nach 
Hilfssatz 1 in zwei positiven Stellen f, < f, geschnitten, und zwar nach (12) unter zwei Winkeln 
%,mit&, >&s. Sei f,im Intervall f, <f; < fs gewählt. Jede Lösung von (5) unter den Anfangs- 
Bedingungen 
Fi)=e, Fi) = as mit %,<S0 So 


kann nach Hilfssatz 1 weder F,(f), noch F,(f) schneiden. Zwei solche Lösungen zu verschiedenen 
&s-Werten würden somit einen Widerspruch zu Hilfssatz 4 liefern. 


7. Schlußbemerkungen 
Der Eindeutigkeits-Beweis gilt auch für $ = 0 ; der Existenz-Beweis aus [1] nicht. Tat- 
sächlich ist ja in diesem Falle nur für © > — 1,2385... eine Lösung vorhanden (vgl. [2]). Diese 
ist dann eindeutig. 
Die Bedingung (7) könnte durch die anscheinend mildere 


Ep ee er; (7a) 


ersetzt werden. Denn es gibt überhaupt keine Lösung von (5), (6), die Werte F(f) > 1 annimmt. 
Eine solche Lösung müßte nämlich mindestens ein Maximum > 1 besitzen, für das die erste 


Gl. (18) die unmögliche Beziehung F > 0 liefern würde. — Auch F(f) = 1 und F(f) = 0 für ein 
endliches f ist unmöglich, da durch diese Bedingungen F = 1 als einzige Lösung von (5) festgelegt 
ist. 

Daß die erste Bedingung (7) für die Eindeutigkeit unerläßlich ist, wurde für genügend große 
positive Werte von CO und ß = lin [3] gezeigt. Ganz anders werden die Verhältnisse bei nega- 


1 
tiven Werten von ß. Darüber unterrichtet die Arbeit [4] für das IntervallO> ß > — > 
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Einige Anwendungen der Theorie der Zylinderschale 
Von Walter Wuest in Göttingen 


ie Theorie der Zylinderschale wird auf die Berechnung nichtgeschlossener Teilstücke erweitert für den Fall 
er en einer Koordlngie Chhan adn Spohsungeteri kan Die Theorie wird dann auf die a 
Biegefeder und die Bimetall-Zylinderfeder angewandt. Die genannte Biegefeder hat eine nichtlineare CI . - 
teristik und wird bei Erreichen einer kritischen Durchbiegung instabil. Im überkritischen Bereich bleibt das 
Biegemoment annähernd konstant und diese Bigenschaft wird zum Antrieb mechanischer Werke ausgenutzt. 


Thetheory of cylindrical shells is extended to open divisions of this shells for the case of uniform stress Her 
bution along one direction. Thistheory isthen applied tothe eross-curved leaf spring and the bimetallie cylinder 
spring. The fore-mentioned leaf spring has a nonlinear characleristic and when reaching a critical point becomes 
unstable. In supercritical statethe bending moment remains approximately constant and this quality is used for 
the construction of spring works. 

La theorie des coques cylindriques est &tendue ü des parties non fermees de celles-ci pour le cas d’une distri- 
bution uniforme des tensions dans une direction. Cette thöorie est ensuite appliqu£e au ressort d lame de section 
courbee ei au ressort bimetallique cylindriforme. Le dit ressort & lame a une characteristique nonlin£aire et en 
atteignant un Etat critigue devient instabil. En Etat supercritique le moment est ü peu pres constant et cetie qualite 
est employee pour la construction de mouvements & ressort. 


Teopna Mammmnpmyeckof IMONOCTU PacHpoctpaunerca Ha PacyeT He3aMKHYTEX YacTeh 
Ha TOT cıyyal, Korza pacıpeneneHnne HANPAKeHUA 3ABHCHT TONBKO OT ONHOK KOOPAHHATEI. 
Naree, Teopma ImpuMmenserca K HONPeYHO-BEILyKIOH IIPY:KuUHe MaTuda uU K ÖHMeTaLImyecKkoNi 
NNNMHNPHYyeckoN IIpysknuHe. YIIOMAHyTaA IIpy>KuHa usruda O61aNaeT He1HHeÜHON XapaKTepHCTH- 
Ko, MH, NOCTUTHYB KPHTUYeCKOrO HPOTH6A, CTAHOBUTCA HeycroüyuBof. B cBeEPXKpHuTnyeckof 
OÖNACTU HsrmÖalımuÜu MOMEHT OCTAETCA IIPHÖNMKEHHO HOCTOAHHBIM,; 3TO CBOÄCTBO HCHONB- 
syeTrcH 1A IIPHBONa MeXAHWSMOB. 


Zur Berechnung ringförmig geschlossener Zylinder wird vonBiezeno-Grammel [l] 
ein allgemeines Lösungsverfahren für beliebige Rand- und Oberflächenbelastungen angegeben, 
das freilich in der praktischen Anwendung eine ganz erhebliche Rechenarbeit erfordert. Einfache 
Lösungen existieren im achsensymmetrischen Fall [2], haben aber nur wenig praktische Bedeutung 
im Hinblick auf technische Anwendungen. Eine größere Mannigfaltigkeit von Lösungen ergibt 
sich, wenn man auch offene Zylinder oder Teilstücke von Zylinderschalen in Betracht zieht. 
Ferner kann man neben den Oberflächenkräften auch eingeprägte Spannungen berücksichtigen, 
wie sie beispielsweise durch wechselnde Temperaturen bei bimetallischer Wandung (Aufbau der 
Wand aus zwei Schichten verschiedener Wärmeausdehnung) auftreten können. Man gelangt so 
zu Problemen der Feinwerktechnik und Meßtechnik, die erfolgreich mit Hilfe der Theorie der 
Zylinderschale behandelt werden können. Die vorliegende Untersuchung ist durch ein Problem 
der Uhrentechnik angeregt worden, nämlich die quergewölbte Biegefeder [3], die als Energie- 
speicher mit annähernd konstantem Drehmoment für technische Uhrwerke besonders geeignet 
ist. Die Untersuchung ist dann auf das verwandte Problem der Bimetallfeder ausgedehnt worden. 
Beide Aufgaben können für bestimmte Belastungsfälle und unter speziellen Randbedingungen 
mit Hilfe der Theorie der achsensymmetrisch belasteten Zylinderschale behandelt werden. Eine 
Erweiterung auf nichtsymmetrische Belastungsfälle ist für eine nachfolgende Veröffentlichung 
vorgesehen. 


Theorie der offenen Zylinderschale bei achsensymmetrischer 
Spannungsverteilung 


Die nachfolgenden Ausführungen lehnen sich eng an die Behandlung der biegesteifen Kreis- 
zylinderschale durch Biezeno-Grammel[1]an. Die Differentialgleichungen der Zylinder- 
schale werden hier nur soweit entwickelt, als sie von den dortigen abweichen und für das Folgende 
notwendig sind. 

Der Hauptunterschied gegenüber der geschlossenen Zylinderschale besteht darin, daß der 
Krümmungsradius «a jetzt nicht mehr konstant sondern belastungsabhängig ist. Ferner sollen in 
die folgende Darstellung auch solche Fälle einbezogen werden, bei denen die Schale bereits im 
spannungslosen Ausgangszustand geringfügig von der Zylinderform abweicht. Bild 1 gibt einen 
Überblick über die benutzten Bezeichnungen, die gewählten Koordinatenrichtungen r, @, z und 
die Verschiebungskomponenten u, ®, w. Mit kyp Kr... UNd Myp My ... bezeichnen wir die 
über die Wanddicke integrierten Spannungen bzw. die Momente dieser Spannungen. Ug, Up %p 
bezieht sich auf die Mittelfläche. 

Die Gleichgewichtsbedingungen gehen im vorliegenden Fall über in: 


dk,, 
ka = a". 2 a Pe: 


apa TR user Biel 


a 
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Wegen der Unabhängigkeit des Spannungszustandes von 9 muß k,, = const sein. Wir setzen 
kz =0, da k, #0 nur einem trivialen Belastungsfall entsprechen würde. Dann ist aber die 
Tangentialkraft %,, gegeben durch: 


k ld 
pp a ( ) 
K 
3 res 
Rage 2322 dz 
Kzz n 
Se Kazt 422 de 
Kort@A2ldz 
dp 
Kope 2K23 do 
Bild 1a Bild 1b 


Bild 1. a Bezeichnung der Koordinatenrichtungen und Verschiebungskomponenten; b Kräfte an einem Ringelement 


Die Momente m,, und m,, stehen in folgender Beziehung zu den Formänderungen der Mittel- 


fläche: 
Br ” Eh? | @ Up | 1 1 ) 4 
Moz B— HEITE, fi u) u da? C,(2) = a m MIay.tıs BI BE ie ( ) 
a Eh?° d2u, 1 1 
MI 75 ae — 69(2) — u 4.) Er 


Dabei ist u = 0,3 die Querkontraktionszahl, c,(2) eine vorgegebene schwache Querkrümmung 
im Ausgangszustand, a, der Längskrümmungsradius im Ausgangszustand und a der Längskrüm- 
mungsradius bei Belastung. 

Aus den Bedingungen (1) und (2) erhält man durch Beseitigen von k,, 


d?m k 
iD BARS Haniy aan me" _ (6) 
und nach Einführen der Beziehungen (3) bis (5): 
d’u 1 d2c 
Sal) le en Re a 


Die quergewölbte Biegefeder 


Die quergewölbte Biegefeder (Bild 2) unterscheidet sich von der normalen Blattfeder nur 
dadurch, daß sie im unbelasteten Ausgangszustand eine Querwölbung aufweist, die mit einem 


Bild 2. Ausgangszustand eines Elementes der quergewölbten Biegefeder 


a 


Bild 3. Gleichsinnige Krümmung durch Biegemomente (Längswölbung Bild 4. Gegensinnige Krümmung durch Biegemomente (Längswölbung 
und Querwölbung im gleichen Sinne) und Querwölbung im entgegengesetzten Sinne) 


Ziehstein in die Feder eingearbeitet wird [4]. Wenn man die quergewölbte Blattfeder entspre- 
chend Bild 3 durch ein Moment belastet, wächst der Biegewinkel nicht linear mit dem Biege- 
moment sondern angenähert parabolisch, d. h. die Biegesteifigkeit des Federquerschnittes nimmt 


FR u 
” rn‘ 
Ei 
= - 
u . 
! 
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mit zunehmender Belastung ab. Bei einer kritischen Belastung, die etwa dem Parabelscheitel zu- 
geordnet ist, kippt die Blattfeder in eine neue Gleichgewichtslage, bei der die Querwölbung fast 
ganz verschwunden ist. Das Biegemoment sinkt dabei auf einen Bruchteil des Anfangswertes und 
wächst bei weiterer Durchbiegung nur ganz langsam an. Man nutzt diese Erscheinung schon seit 
langem bei Bandmaßen aus, die in gestrecktem Zustand sehr biegesteif sind, sich aber im ge- 
krümmten Zustand leicht zusammenrollen lassen. In neuerer Zeit benutzt man die querge- 
krümmte Biegefeder wegen des konstanten abgegebenen Momentes (im überkritischen Fall) zum 
Antrieb technischer Uhrwerke. m. 

Im folgenden wird gezeigt werden, daß wesentliche Eigenschaften dieser Feder bereits mit 
der einfachen Annahme achsensymmetrischer Spannungen berechnet werden können!). Unter der 
Voraussetzung einer schwachen Querwölbung (praktisch erhält man noch eine brauchbare Nähe- 
rung, wenn b/r, nicht wesentlich größer als 1 ist) und einer Beanspruchung durch ein reines Biege- 
moment kann die deformierte Blattfeder als Teil einer schwach gestörten Zylinderschale ange- 
sehen werden, wobei die Achse des Zylinders in Richtung der schmalen Seiten in Bild 3 bzw. 4 
verläuft. Wenn die Spannungsverteilung achsensymmetrisch, d.h. unabhängig von 9 sein soll, 
müssen die Biegemomente in ganz bestimmter Weise an den schmalen Seiten angreifen. Dies 
wird im allgemeinen nicht der Fall sein, doch kann man unter der Voraussetzung, daß die Länge 1 
der Blattfeder groß gegen ihre Breite 5 ist, nach dem Prinzip von St. Venant annehmen, daß die 
hierdurch hervorgerufenen Störungen auf die unmittelbare Umgebung der Lastangriffsstellen 
beschränkt bleiben, während der überwiegende Teil der Blattfeder einer achsensymmetrischen 
Spannungsverteilung unterworfen ist. Von den genannten Einschränkungen abgesehen, beschrei- 
ben die oben aufgestellten Gleichungen den Spannungs- und Deformationszustand der Blattfeder. 
Dabei ist in Gl. (4) und (5) c, = 1/r, und a, = ® zu setzen, wenn vorausgesetzt wird, daß die 
Feder im Ausgangszustand gestreckt ist und eine kreisförmige Querkrümmung mit dem Radius r, 
besitzt (Bild 2). Gl. (7) geht also über in: 

da 
am 0 en a 
mit 
a? 
nl —A) 72- 


Wenn man gleich berücksichtigt, daß u, eine symmetrische Funktion in z sein soll, lautet 
die Lösung: 


z 2 a z 2 
GC Era Ten > Nena —- san nee re : 
= C,Cojn — cosn +0, &inn “sinn (D 


Die Integrationskonstanten (©, und 0, bestimmen wir durch die Bedingungen: 


& b S b 

kr =0 für :=4+7Z: MM =0 für Be dee ERELIDERER EL 
Durch Ausrechnen dieser Bedingungen folgt für die Integrationskonstanten: 
nb . nb .%b nb 
ee ea 

h rn HK a)Im De e nbzienbunpen Ser Perl: 
2a 3. rer 
nb , nb 9b nb 
lh (d 24) Mr Co, inz, + Sin z— cos. 

2 rg ar en m nb . nb ee ee 

Ta TE ya 


Für die praktische Anwendung interessiert vor allem die Abhängigkeit des Längsbiege- 
momentes M von der Längskrümmung. Das Gesamtmoment M ist aber gegeben durch: 


B/2 
M =) Moe + Konto] de TR REN 


1) Anmerkung bei der Korrektur: Eine Näherungstheorie der quergewölbten Biegefeder i i 

. +1 1» * t b t 
P. Funk: Über ein Stabilitätsproblem bei den durch Krümmung steif eniesigen Meßbändern. Osten 
Arch. V (1951), H.4,S. 387—397 gegeben worden. Dabei ist allerdings die Annahme zugrunde gelegt worden, 
daß der Querschnitt kreisbogenförmig bleibt und sich nur der Querkrümmungsradius bei Belastung ändert. ; 


And u 


Z. angew. Math. Mech. N : B 
Bd.34 Nr. 12 Dez. 1954 Wuest, Einige Anwendungen der Theorie der Zylinderschale 447 
ee EEE nen ei RT Eu ae eure ee ae Eric 


Die Auswertung dieses Integrals, die im Anhang wiedergegeben ist, führt zu folgendem Ergebnis: 
3 
M ze +u2)—u (* +) Pi+(E +) P, Br, 
Die Funktionen ?, und P, sind dabei gegeben durch: 
_ 2 @ojA — cosA 
ent EsinA 
p2 1 Cofja—cosi _ Sin Asin 
24 SinA+sini (SinA+sinA) 


nb a4 ur En oe 
== — ya — 
A 3( Sr 


Hierbei ist 


a 


Nee P, und P, sind in Tabelle 1 berechnet und in Bild 5 aufgetragen. Wir führen 
noch ein: 


Bild5. Hilfsfunktionen ?P, und P, zur Berechnung Bild 6. Hilfsfunktion F zur Berechnung des Biegemomentes 
des Biegemomentes nach Formel (14) nach Formel (15) 


wobei also jetzt A, für ein bestimmtes Federband eine Konstante ist. Wir können dann für das 
Moment schreiben: 


_ bA® ZAS >) 
M ee 5. (1 — u) F (% F} a . . . . . . . . . . . (15). 
Dabei ist F nach (14): 
a Rz) ir rau 


Die Funktion F ist in Bild 6 für verschiedene Werte von A, aufgetragen. Man erkennt, daß diese 
Funktion und damit auch das Biegemoment einen Größtwert erreicht, der dem kritischen Fall des 
Abkippens entspricht. 


Tabellel. Funktionen P,() und P,(A) 


%» | Pia) Bü) | %- | „Pi P,(A) 
0 1,00000 0,00000 3,0 0,72564 0,16771 
0,5 0,99968 0,00035 3,5 0,61792 0,17683 
1,0 0,99450 0,00538 4,0 0,52692 0,16108 
1,5 0,97296 0,02600 5,0 0,40372 0,11420 
2,0 0,92112 0,07001 6,28 0,3174 0,0793 
2,5 0,83434 0,12768 
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Dieser Maximalwert und damit auch die kritische Krümmung ist offenbar nur von Ay ab- 
hängig. In Bild 7 sind die kritischen Werte r,/a, und F in Abhängigkeit von A, aufgetragen. 


Tofax Fk 


Vergleich mit Messungen 
Für die Messungen stand ein Feder- 
band mit folgenden Abmessungen zur Ver- 
fügung: 
b=20mm h=02mm n=16mm. 


Daraus folgt 
a 14,35 
ki 12 Va = 14,35. 
Diesem entspricht nach Bild 6 


0 — 0,032 
477 


und damit 
Bild 7. Zur Berechnung des kritischen Längskrümmungshalbmessers a7, 


und des kritischen Biegemomentes a; —=500 mm. 


Gemessen wurde in guter Übereinstimmung damit a; = 490 mm. Die bisherigen Ergeb- 
nisse können näherungsweise auch auf den Fall übertragen werden, daß die Blattfeder an einem 
Ende fest eingespannt ist und am anderen Ende eine Kraft wirkt, so daß also das Biegemoment 
und auch die Steifigkeit längs der Blattfeder veränderlich ist. Da auch der Längskrümmungs- 
radius a längs der Feder veränderlich ist, kann man diese nur abschnittweise als Zylinderfläche 
auffassen. Auch für diesen Fall ist eine befriedigende Übereinstimmung zwischen den gemessenen 
und den näherungsweise berechneten Werten gefunden worden (vgl. [3]). Doch ist dabei wesent- 
lich, daß das feste Ende der quergewölbten Blattfeder nicht zwischen ebene sondern gewölbte 
Flächen eingespannt wird. Solche Fälle dagegen, bei denen die Einspannung von wesentlichem 
Einfluß ist (z. B. kurze Federn) können nicht mit der achsensymmetrischen Theorie behandelt 
werden. 


Überkritischer Zustand 


Bei Überschreiten der kritischen Größen kippt die Feder in den überkritischen Zustand, bei 
dem die Querkrümmung bis auf eine schmale Randzone ganz verschwunden ist. Im Fall der 
gleichsinnigen Krümmung (Bild 3) erreicht die Feder schließlich ebenfalls den überkritischen Zu- 
stand; doch geht dieser Übergang allmählich und ohne Kipperscheinungen, dafür aber unter Tor- 
sionserscheinungen, vor sich, die rechnerisch nicht leicht zu erfassen sind. 

Man übersieht leicht, daß die Gl. (1) bis (5) auch die triviale Lösung u, = k,, = kyp = 0 
haben, woraus folgt: 


RR R] 1 1 

Mg = 775 E +u = aale id) 
BE 1 1 a) 

Mat iD —G +u rer (17). 


Diese Lösung erfüllt aber nicht die Randbedingung M;,, = 0 fürz = + b/2 und gilt also zumindest 
nicht in der Randzone. Wenn man dies nicht beachtet, folgt für das Gesamtmoment M anstelle 
von (14): 


M = 


BIEBE: (2 s) solltet 


Be rn a 1a To 

Formel (18) ist aber nichts anderes als der Grenzwert von Gl. (14) für große Werte von). 
Formel(18) gilt übrigens auch für den Fall gleichsinniger Krümmung (Bild 3), wenn man nur den 
Krümmungsradius «a negativ einsetzt. 

Formel (18) kann dadurch nachgeprüft werden, daß man die quergewölbte Blattfeder ent- 
sprechend Bild 8 belastet. Wenn man die geraden Schenkel lang genug wählt, wird der halbkreis- 
förmige Abschnitt durch ein annähernd konstantes Moment belastet. Es ist dabei auffallend, daß 
der sich einstellende Krümmungsradius gerade gleich dem anfänglichen Querkrümmungsradius r, 
ist, Dies kann man auch leicht rechnerisch nachweisen, indem man für einen zunächst beliebig 


£ 


i 
: 
E 


Danger, Math. Mech __ Wüuest, Einige Anwendungen der Theorie der Zylinderschale 449 


angenommenen Krümmungsradius a die aufzuwendende Formänderungsarbeit berechnet2). Da es 
sich im vorliegenden Fall um einen reinen ebenen Hauptspannungszustand handelt, ist die be- 
zogene Formänderungsarbeit gegeben durch: 


G p 
Ar —— 1m (e? u & — 2 uE1&) . 
Die gesamte Formänderungsarbeit erhält man durch Integration über 2 
die Blechdicke, die Breite und einen Halbkreis zu 
G bh’n | | 1 1 1 
A=— Be an lie 
De ae lea er, . 


. . . D Re .. Bild 8. Messung der Federsteifigkeit 
Der Ausdruck in der eckigen Klammer ist aber ein Minimum für rat a Komet 


a—-+r, d.h.im Fall vonBild8ist der Längskrümmungsradius gleich (18 mit «= +,). Die Messung 
d fä lich krü di daß di wird einmal mit gleichsinniger und 
em anfänglichen Querkrümmungsradius, vorausgesetzt, daß die ge- einmalmit gegensinniger Krümmung 


raden Schenkel lang genug sind, so daß im gekrümmten Teil nahezu eeanz 
ein reines Biegemoment übertragen wird. 


Zur Nachprüfung von Formel (18) sind zwei Federn entsprechend Bild 8 belastet worden. 
Im Fall der gleichsinnigen Krümmung ist ein geringeres Gewicht erforderlich, um den oberen 
Schenkel in der Schwebe zu halten als im Fall der gegensinnigen Krümmung. 


Beder Ab 16mm 3% = 0,10 mm; +, = 16 mm- 
Setzt man ferner E = 2100 kg/mm?; u = 0,3; a =r, so wird das Biegemoment: 
|M| = 1,92 (1+ 0,3) kg mm 


wobei das Minuszeichen sich auf den Fall der gleichsinnigen Krümmung bezieht. Gemessen wurde 
demgegenüber: 


|M| = 2,16 (1 + 0,27) kgmm 
wenn für M der Mittelwert des Momentes im halbkreisförmigen Abschnitt genommen wird. 
Feder B: b = 20,4mm; h = 0,20 mm; r, = 16mm 
|Mineor.| = 15,40 (1 + 0,30) kg mm 
IM gemess.| = 14,35 (1 + 0,29)kgmm. 
Die Übereinstimmung ist also in beiden Fällen befriedigend. 


Beanspruchung 


Die maximalen Biegespannungen sind im überkritischen Fall, abgesehen von der Rand- 
zone, gegeben durch: 


E ı/l 4 
4 | > TO. hen Bsissiuft Anassenm: ah 
Pa 1— u? tn, (19) 
E 4 y 
—= -— en ee 
9 1— u? 2\a pi No a 
In der Randzone fällt dagegen o, auf den Wert O0 ab, während 
E A ee 1 En Up 
en u 


ist. 
Nun ist aber nach (9), (11) und (12): 


se ( An} (2%) wie 
“(7 ee a) en ee 


so daß die tangentiale Spannung am Rande für große Werte von A gegeben ist durch: 


ee Babe Bean Ya a 
21a Y3 (1 — 2) To a 


2) Anmerkung bei der Korrektur: Vgl. hierzu auch: P. Funk: Über die durch Krümmung steifgemach- 
ten Meßbänder. Z. angew. Math. Mech. 14 (1934), H.4, S. 251—252. 


ER A 
| | i % 


b+20 he02 1* mm 6 
I 95cm lang gegensinnig | 
IT 35 cm lang gleichsimg | 
II 35cm lang Differenz I-I 
IV 347cm lang gegensinnig 


a 
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Bild 9. Messungen an spiralförmig in einem Federhaus aufgewickelten quergewölbten Bandfedern 


allen Fällen steigt das Drehmoment nur wenig mit dem Drehwinkel an. Bildet man die Differenz 
zwischen dem Fall gegensinniger Krümmung (I) und gleichsinniger Krümmung (II), so ist die 
Differenz (III) nahezu konstant. Das Differenzmoment ist aber leicht aus Formel (18) zu be- 
rechnen und offenbar gegeben durch: Aut 


a ne re. 


Für das gemessene Beispiel (Bild 9) findet man M) = 1,16 kg cm, was nur wenig unter dem ge- 
messenen Differenzmoment von 1,2—1,3 kg cm liegt. Wenn man die sonstigen gemessenen Kur- 
ven mit Formel (18) vergleicht, folgert man, daß im praktisch ausnutzbaren Bereich r,/a zwischen 
0,8 und 1,3 schwankt. Dabei wird bei der kurzen Feder der Schwankungsbereich bereits nach 
wenigen Umdrehungen durchlaufen, bei der langen Feder erst nach einer entsprechend größeren 
Anzahl Umdrehungen. 


Ebene Blattfeder 


Wenn man eine im Ausgangszustand völlig ebene Blattfeder durchbiegt, zeigt sich eine 
schwache Querwölbung, die durch die Querkontraktion hervorgerufen wird. Man kann diesen 
Effekt besonders deutlich an einem Gummistreifen beobachten, weil dort die Querkontraktion 
groß ist. Dieser Querwölbungseffekt ist in den vorher abgeleiteten Gleichungen mit enthalten, 
wenn man nur 7, = setzt. Die Formel (14) für das Biegemoment geht damit über in: 

bP® ‚1—wWP,+tuP, DbRM 


hit Te 12a 


Die Funktion G,(A) ist für u = 0,3 in Bild 10 aufgetragen. Dabei ist 
ne al: 
= )3(1— u2) | — 
ı-Ba—e 7, 
b Breite der Blattfeder, } Blechdicke, a Biegeradius. 
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Wir wollen ferner noch den Querwölbungsradius r„, in der Mitte der Blattfeder bestimmen 
den man durch zweimaliges Differentiieren von Gl. (10) an der Stelle 2 = 0 findet. Es ergibt sich: 


Co & sin + Sin, cos 


Im TEENS AT ET = u G, (A) AST TER FRE 


Die Funktion @,(A) ist ebenfalls in Bild 10 
aufgetragen. Man erkennt daraus, daß die 
Querwölbung nur bei kleinen A-Werten 
sichtbar sein wird. 


Bild 10. Hilfsfunktionen @, und G, bei der ebenen Blattfeder nach Bild 11. Wärmedehnungen und elastische Dehnungen in 
Formel (23) und (24) einem Element des Bimetallbleches 


Bimetallblech 


Ein Bimetallblech besteht aus zwei Schichten von Metallen mit verschiedener Wärmeaus- 
dehnung, so daß sich ein Bimetallstreifen bei Wärmeeinwirkung durchkrümmt. Obwohl an der 
Grenzfläche der beiden Metallschichten durch den Herstellungsgang ein Verschweißen eintritt, 
wollen wir doch diese Grenze als Unstetigkeitsfläche der Materialeigenschaften idealisieren und 
die beiden Metallschichten als homogen voraussetzen. 

Da die Spannung o, in Richtung der Wanddicke entsprechend den Annahmen der Schalen- 
theorie verschwindet, folgt für die anderen beiden ee 2 


(tu); en (&-+us)- 


O4 = 
1 Zi 


Far 


Bezeichnen wir mit &,, und &,, die Dehnungen an der Grenzfläche und mit x, und x, die Krüm- 
mungen, so sind die Dehnungen in den beiden Richtungen: 


g=E&o +21 —aT, S=tpt 12 —üT. 


Dabei ist x an der Grenzfläche gleich 0 und läuft von —h;; bis +h, (vgl. Bild 11). Für a, ist je 
nach der Schicht, in der man sich befindet , oder a,; einzusetzen. T ist dabei die Temperatur- 
änderung gegenüber der spannungsfreien Ausgangstemperatur. 


Durch Integration über die Wanddicke ergibt sich für die Kräfte und Momente: 


Exhr+ Ey (nen) + (Brt— Bucht) (2 +0) 
= 0 See nn» (28), 


(Erhı — Errhrr) (eo + MEor) + (Erhit+ Er kn) z 2 (m + 4%) 
— (BE, hir — Enhtro) 1 + er T=—(l—u)2m;s - - - (26)- 


j,k= 1,2 bzw. 2,1 entsprechend der Längs- und Querrichtung. (25)und (26)'stellen vier Gleichungen 
zur Bestimmung von gg}, &gg, %ı und x, dar. 
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Ebenes Bimetallblech mit geringer Auslenkung vum 
Wenn das Bimetallblech keinen äußeren Kräften und Momenten ausgesetzt ist, gilt in diesem 
Fall k,., = ka = Mia = My = 0. Damit erhält man: 


Ei = ©0253 %ı = Ka . 


Dies führt aber zu den vom Verf. in zwei früheren Veröffentlichungen [5] abgeleiteten Beziehungen 
für den Bimetallstreifen, auf die hier nicht näher eingegangen werden soll. 


Bimetall-Zylinderblech 
Zylinderförmig gebogene Bimetallbleche werden vor allem in meteorologischen Registrier- 
geräten angewandt. Dabei hat es sich gezeigt, daß breite Zylinderbleche einen merklich größeren 
Ausschlagwinkel bei Temperaturänderungen zeigen als schmale Streifen aus dem gleichen Mate- 
rial [6]. Dieser Effekt soll im folgenden berechnet werden. Dabei beschränken wir uns auf den 
Sonderfall Z=E,, =E und hr =Ähr =hj/2, der bei den Anwendungen in der Meßtechnik 
meist angestrebt wird. Mit dieser Vereinfachung geht die erste Gl. (26) über in: 


12 3 ger’, 
(1) 75 Me ta) Zum) ATMT- 


Wir wenden auch auf dieses Problem die Theorie der symmetrisch belasteten Zylinderschale an. 
Dabei ist allerdings zu bedenken, daß die Voraussetzung der Näherung, nämlich abgewickelte 
Federlänge groß gegen Federbreite, in praktischen Fällen viel schlechter erfüllt ist als bei der 
quergewölbten Bandfeder. Davon abgesehen sind die Gl. (1)—(7), auch hier gültig, wenn 
Co (2) = 0 gesetzt wird und die Indices 1,2 mit o, z gleichgesetzt werden. Die weitere Rechnung 
geht in gleicher Weise vor sich wie bei der quergewölbten Biegefeder und man erhält als Ergebnis 
statt der Gl. (14) für das Gesamtmoment: 


Be ee .- T r+% 
E 5 IM FR Hz Ted) — ".G-,) 
3 T 1 1 3 Tr]? 
nt re) Feet er... en. 


Dabei sind P, und P, die gleichen in Bild 5 aufgetragenen Funktionen. Setzt man M =(, so 
kann man den Ausschlag eines nichtbelasteten Bimetallbleches unter der Einwirkung von Tempe- 
raturänderungen berechnen. Bei Beschränkung auf kleine Temperaturänderungen kann das 
quadratische Glied in (27) vernachlässigt werden und man erhält für die Änderung Aa des anfäng- 
lichen Halbmessers ag: 


da _ 3 a „d+mW)d—uB)) 
U — 2 (a a &ır) ı AR 37 Tempe = Kun a. es Fre Se Fleiete (28). 
7 Die Funktion P,(A) ist in Tabelle 1 tabelliert. 
Für A>5 gilt näherungsweise P,(A) = 2. 
13 Dabei ist 
4 2 2 
= BA—m)YLr1285 1/2. 
12 ha ha 
Der Vergrößerungsfaktor 
1 r,_+W4-—uP) 
Lee Bu 
ı um den sich der Ausschlag eines breiten 
e 5 10 15 20 25 Bimetall-Zylinderbleches gegenüber einem sehr 
en . schmalen vermehrt, ist in Bild 12 aufgetragen. 
Biete Drttene BD ” Der Ausschlag kann bis zu 30%, größer werden. 


Zusammenfassung 


Während die Theorie der symmetrisch belasteten Kreiszylinderschale nur wenig Anwen- 
dungen gefunden hat, solange man sich auf ringförmig geschlossene Schalen beschränkt, gibt es 
verschiedene technische Probleme, die mit der gleichen Theorie, jedoch als offene Schalen behan- 
delt werden können. Hierzu zählen die quergewölbte Biegefeder, die als Kraftspeicher von an- 
nähernd konstantem Drehmoment in der Uhrenindustrie Verwendung findet, sowie das Bimetall- 


N 
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Zylinderblech, das als Meßelement in Temperatur-Registriergeräten dient. Für beide Probleme 


sowie für den Wölbungseffekt bei der gewöhnlichen Biegefeder werden Berechnungsgrundlagen 
hergeleitet. 


Anhang 
Berechnung der Funktionen P,(A) und P,(A). 
Durch Einsetzen von (3), (4), (6) und (7) in (13) folgt mit a, = &: 
bj2 
eb na 1 1 = | h \ n 
a a al ba 
Demnach ist 


1 1 wer = 
ee 
y 1 


4 Ey) A } 
— 2r fc, + 0) Ein & cos 4 — (0, — 0) Co sin 2] 


Sin a cos? = + 60? -- sin? Z 


. z Co A Sin A 5 + eos = sin u 
2 2 
mit A 2 s 
Daraus folgt: 
2 &0jA — cos} 


' A 6©ini+sini' 
Entsprechend ist: 
i b/2 


bh 1 we Yp [" 
— [ud 
ra zen ACH nur Sr nie Me 
—b/2 
b/2 
Nn2 


ß Cop"? — cos? — no 2.0. Co - Sin 2? cos — sin 7 
5/2 


I 


NZ 


+60: Sin"? sin? | dz 


| 


| & ua (Cor er Bea \ [ G— 0 Cof Be cos Be 1 
a 


2 2 2 
9/2 
C,6, ANZ und | 
+ — 9 Sin — 5 sin —- de. 
Nun ist aber 
b/2 9 
(607 — sin? nn) de == (Sina — sinA) 
j a a a 
—/2 
v2 } 
| (601° cos +1) 1) de = 2" (Sina c0s2 + Cojasina) +D 
[2 _ 
bj2 


zZ 2anz 


: Ein - sin — dd= 5 CojA sind — Sin cosA) 


—b/2 
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10 = bie ah = CojA — cosA 


7 a) \2n2) “(SinA + sinA)® 
1 1\2/ a2\2 SinA sin} 
Tr ai“ ) (6) (Ein + sinA® 


2m? 


1 1\2/ a\2. 1—EojA cosi 
010 --(, u =) | (Sin + sinA)' a 
Zusammenfassend ergibt sich 
/ 
ie (4 Ay a Gaumen N Sinasind 
[“ 2 ar 21&inA+ sinA (SinA + sinA)2l 
= “ 


Daraus folgt aber 


8n? 


1 &ojA — cos} SinA sinA 


Ten ne 
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Ein neues Auswerteverfahren der ebenen Spannungsoptik 
Von L. Föppl in München 


Der ebene Spannungszustand wird auf das orthogonale Netz der Hauptspannungslinien bezogen. Die 
Gleichgewichtsgleichungen werden mit Hilfe der zur Spannungssumme als Realteil gehörenden komplexen 
analytischen Funktion dargestellt und integriert, wobei der Integrationsweg in der Ebene beliebig gewählt werden 
kann. Der Integrand enthält dabei die auf spannungsoptischem Weg an jeder Stelle gewonnenen Werte der . 
Differenz der beiden Hauptspannungen und ihre Richtungen. 


The orthogonal lines of principal stress are used as co-ordinate lines to describe the plane state of stress. T’he 
conditions of equilibrium are represenied byaid of an analytic function that has the sum of the principal stresses 
as real part, and they are integrated along an arbitrary path in the complex plane. The integrand depends on 
the difference of the two principal stresses and their directions, that are found by photoelastic means in any point. 


L’&tat plat de tension est refere au reseau orthogonal deslignes de tension principales. Les &quations d’equilibre 
sont presentees et integrees a l’aide de la fonction complexe, analytique, correspondant & la somme de tension 
comme partie reelle, ü quoi le chemin d’int&gration dans la plaine peut ötre choisi arbitrairement. A cela 
l’integrand contient les valeurs de la difference des deux tensions principales, obtenues par la voie de tension 
optique a chaque endroit, et leur directions. 


IlsockocTHoe HANPAKEHHOE COCTOAHUe OTHOCHTCHA K CETU TJIABHEX TPaeKTopuH Hanpa- 
;keHuus. PaBHoBecHsie YpaBHeHnAa IPeNCTaBıAmTcH IIPH IIOMOINM KOMILIEKCHEIX AHANHTHYSCKUX 
PyHkunA, TPUHAAISKAUMX, KAK BEIIECTBEHHAA YACTb, K CYMMe HANPA:KeHHÜ, U HHTETPHPYIOTCA 
Ipu IPOHSBONbBHO BbIOPaAHHOM KOHTYpe B I1OCKocTn. [lonuHrerpansHoe BErpaskeHue CoNepsKHT 
3HayeHHA PA3HOCTU O6OUX TJIABHEIX HAIPFIKEHHÄ U UX HANPABIeHHÄ, IOMYYeHHLBIe B KAKoli 
TOYK® H3 OIITHYECKUX SbheKToB HAanpsskKeHng. | 


Einleitung 

In der ebenen Spannungsoptik erhält man durch eine erste Aufnahme das Bild der Isochro- 
maten, längs deren die Differenz der Hauptspannungen o, und o, konstant ist. Nachdem Polari- 
sator und Analysator samt den damit verbundenen A/4-Plättchen miteinander vertauscht worden 
sind, wodurch aus einem Circular-Polariskop ein Linear-Polariskop entsteht, erhält man durch 
eine zweite Aufnahme die zu dieser Stellung des Linear-Polariskopes gehörige Isokline. Sie 
verdankt ihre Entstehung der dem Linear-Polariskop anhaftenden ausgezeichneten Richtung, 
indem die am Modell ankommenden Strahlen in einer bestimmten Richtung linear polarisiert 
sind. Alle Punkte des ebenen Spannungszustandes, in denen eine der beiden Tangenten an die 
Hauptspannungslinien mit dieser ausgezeichneten Richtung übereinstimmten, lassen den in 
dieser Richtung polarisierten Strahl durch, ohne daß sich in der dazu senkrechten Richtung, 


ie un ' um de 
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die zur anderen Hauptspannungslinie in diesem Punkt des ebenen Spannungszustandes gehört, 
eine Teilschwingung ausbildet. Nach dem Durchgang durchs Modell bleibt demnach der Licht- 
strahl in dieser Richtung linear polarisiert und wird vom Analysator, der nur eine senkrecht hierzu 
schwingende Lichtkomponente durchlassen würde, vollkommen verschluckt. Die Verbindungs- 
linie aller dieser Punkte nennt man eine Isokline, die als dunkle Linie neben den als Jsochro- 
maten gekennzeichneten dunklen Linien auftritt. Im Gegensatz zu letzteren wandert die Jso- 
kline bei gleichmäßiger Drehung der beiden Polarisatoren der spannungsoptischen Apparatur 
über das Gesichtsfeld. Indem man jeweils etwa nach 10° Drehung eine Aufnahme macht, erhält 
man die zugehörige Isokline und indem man all’ diese Isoklinenaufnahmen auf ein Zeichenblatt 
zusammenstellt, erhält man das Feld der Isoklinen. Längs jeder Jsokline hat das 
Achsenkreuz der Hauptspannungslinien im Raum die gleiche Richtung, daher der Name ,Jsok- 
line“. Dem Übergang von einer Isokline zur nächsten entspricht jeweils eine bestimmte Drehung 
dieses Achsenkreuzes um einen Betrag von 10°, wie wir oben beispielsweise angenommen hatten. 
Indem man diese Achsenkreuze längs aller Isoklinen einträgt und sie zu Linien verbindet, erhält 
man das Orthogonalsystem der Hauptspannungstrajektorien. Mit der Jsochromatenaufnahme 
und dem aus dem Feld der Isoklinen gewonnenen Netz der Hauptspannungslinien erschöpft sich 
die einfache ebene Spannungsoptik. Sie liefert demnach nur die Richtung der beiden Hauptspan- 
nungen o, und o, an jeder Stelle und deren Differenz o, — o,. Dagegen kennt man damit noch 
nicht die Größe der Hauptspannungen o, und o, einzeln. 

Es gibt verschiedene Verfahren, um aus der Kenntnis der Differenz der beiden Haupt- 
spannungen an jeder Stelle sowie der Kenntnis des Netzes der Hauptspannungslinien die Haupt- 
spannungen selbst zu bekommen. Es hat sich herausgestellt, daß keines der bisher bekannten 
Verfahren zur vollständigen Auswertung des ebenen Spannungszustandes immer brauchbar wäre. 
Es ist deshalb von praktischer Bedeutung, wenn man noch weitere Auswerte-Verfahren ent- 
wickelt. Im Folgenden soll ein neues Auswerte-Verfahren angegeben werden. 


$ 1. Die Gleichgewichtsbedingungen und die Verträglichkeitsbedingung 
In Bild 1 sind die gegen das raumfeste x-y-Koordinatensystem um den Winkel «x bzw. 
a 25 gedrehten Linienelemente ds, und ds, der Hauptspannungslinien eingetragen. Der 


Neigungswinkel a der Hauptachse 1 gegen die x-Achse ist eine Funktion des Ortes. Die Kurven 
& = constin der x-y-Ebene, in der sich der zu untersuchende ebene Spannungszustand befinden 
soll, entsprechen den Jsoklinen. Aus dem in der Einleitung besprochenen Feld der Jsoklinen 
läßt sich a(z, y) an jeder Stelle des ebenen Spannungszustandes entnehmen. 


Bild1 Bild 2 


Wir wollen die Spannungen auf das orthogonale Netz der Hauptspannungslinien beziehen !). 
Dabei gehen wir vom rechtwinkligen Koordinatensystem x, y aus (siehe Bild 2). An einem drei- 
eckigen Scheibenelement, das von zwei Hauptspannungslinien mit den Linienelementen ds, und 
ds, begrenzt wird und dessen dritte Seite parallel der x-Achse verläuft, greifen die Hauptspannun- 
gen 0,, 0, sowie die Normalspannung o, und Schubspannung 7 an. Aus dem Gleichgewicht in 
Richtung der x- und y-Achse folgt: 


oy=0,5in?« + 0,c05?0; 
T = (0, — 05) Sina Cosa; 


1) Siehe z.B. L. Föpplu. H. Neuber „Festigkeitslehre mittels Spannungsoptik“ München, Eerlin 
1935, 8.13. 


456 Föppl, Ein neues Auswerteverfahren der ebenen Spannungsoptik pa.34’Nr.12 Dez. 1054 


und entsprechend gilt 
0, = 0,6089? & +4 0,5in?a. 
Mit den Abkürzungen: 


ae hd BR a ee (la), 
2 
TE man At . (1b), 
lassen sich die obigen Gleichgewichtsbedingungen umschreiben: 
Y=P9 HI:0820) als Br (2a), 
ed LITE COSZU NE W LlE SCHUNEE ee (2b), 
Tu. BD Zt Aue er u ee (2c). 
Setzt man diese Spannungswerte in die Gleichgewichtsgleichungen: 
00T 
N, PRPFRALANEIE 200 ER 3a), 
02 0 x > 
00, 90T 
ae 1 ae Krk 3b 
au T dm (3b) 
ein, so erhält man mit den Abkürzungen: 
hasYy)=g.)0082%, las) gg sn 2% 
op oh 
=... .... (4a), 
0% 02. 0% ei 
op af A, 
EEE LEN 


Bevor wir auf diese Gleichungen zurückgreifen, wollen wir noch die Verträglichkeitsbedingung 
beim vollkommen elastischen Spannungszustand ausnützen. Sie läßt sich auf die Gleichung 


zurückführen. Um p, das die Dimension einer Spannung hat, dimensionslos zu machen, denken 
wir es durch den Elastizitätsmodul Z dividiert. Aus Gründen, die sogleich Erwähnung finden wer- 
den, dividieren wir p durch E/2. Die Funktion 2 p/E ist nämlich der sogen. „Drehung“ 


> ale 2.) ar ak En (6) 


an jeder Stelle des ebenen Spannungszustandes konjugiert zugeordnet, so daß 
= +tio= K(z) 


eine analytische Funktion K(z) der komplexen Veränderlichen 2=x-+y darstellt. Fürs 
Weitere wollen wir den Faktor 2/E von p weglassen und nur in Gedanken hinzufügen, so daß wir 
die letzte Gleichung vereinfacht so schreiben wollen: 


ORT N een (7) 


Was die Definition der elastischen Drehung ® nach Gl. (6) betrifft, so stellt sie die mittlere 
Drehung der beiden Linienelemente dx und dy entgegengesetzt dem Uhrzeigersinn dar. Mit € 
und n sind dabei die elastischen Verschiebungen parallel der &- bzw. y-Achse bezeichnet. 
Die Drehung @ ist eine reine Ortsfunktion, ebenso wie p und genügt auch der Laplaceschen 
Differentialgleichung 

a ED RT Er Re VS (8). 


Da K(z) eine analytische Funktion von z ist, gelten de Cauchy-Riemannschen Glei- 
chungen: 


op _ 
au en ae) 
op [N] 


%y = I EEE re a ae N (9b), 
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oder auf das rechtwinkelige Netz der Hauptspannungslinien mit den Linienelementen ds, und ds, 
bezogen: 


op dw 

ös, = ds, ETF ee (10a), 
0 0m 

Fee npatal ab. vb Samui sa Sl (10b) 


Schließlich sei noch darauf hingewiesen, daß man das Linienelement in der Ga uBß schen Ebene 
auch in Polarkoordinaten umschreiben kann: 


d&=de-+idy=ds-cosp-tids-sinp=ds-e? ...... (11); 
wenn mit p der im Gegenuhrzeigersinn positiv gezählte Winkel bezeichnet wird, den das Linien- 
element ds mit der positiven x-Achse einschließt. 

Wir gehen von Gl. (7) aus und differentiieren die komplexe analytische Funktion K(z) 

nach z: 

er : ‚io 0p .0p 

een a a Fr 

Wir haben dabei Gl. (9b) berücksichtigt und davon Gebrauch gemacht, daß der Differential- 

quotient einer komplexen analytischen Funktion von der Richtung, in der differentiiert wird, 

unabhängig ist. K,(z) ist wieder eine komplexe analytische Funktion. Unter Beachtung der 
Gleichgewichtsgleichungen (4a) und (4b) kann man sie umschreiben: 


u. 

K,() = (e + öy —1i Tun ELCH A Are (13) 
oder, indem man die Abkürzungen f, und f, eliminiert, 

N 0q =): 

KBo=— Al + 20%) 42-205, 


Es sei noch darauf hingewiesen, daß die Gl. (4) auf das rechtwinkelige Netz der Haupt- 
spannungslinien mit den Linienelementen ds, und ds, bezogen werden können, indem man das 
x-y-Koordinatensystem um den Winkel « gedreht denkt, so daß die x-Achse in die Rich- 
tung des Linienelementes ds, und die y-Achse in die von ds, fällt; damit folgt aus den Gl. (4) 
mit de—ds,; dy—ds, unda—=0: 


MINI IR V TRBDIERURSEN WEHR 192), 
Os, 05} 085 ee 
op _0q En (15b) 
05 08 %s 
oder 
00] 0% 
], ZEIT DIR TAU AQ3S) 1 DE DR 16a), 
8, I, 
005 0% 
—e ter a er BR ee} AR (166), 
ds, 2gq 35, (16b) 


Wenn man noch die Krümmungsradien R, und R, der beiden Hauptspannungslinien am be- 
treffenden Punkt einführt: 


= E er er reise u; (17a), 
= ® DA RE PEN PR (17b), 
so erhält man die bekannten Gleichgewichtsbedingungen 
= ar ® Pan» ve FE Faro BP AN RERER (18a), 
we = ” ne A ac: (18b) 


— 92 z1s reine Ortsfunktion auftritt, kann Gl. (18b) aus Gl. (18a) nicht durch 


o 

Damit g= -— 

Vertauschung der Indizes 1 und 2 gewonnen werden, weil sonst g sein Vorzeichen wechseln würde, 
sl 


u u au 
4 er. - 
Z. angew. Math. Mech. 
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$ 2. Das Auswerte-Verfahren 

Wir greifen auf Gl. (14) zurück und integrieren sie zwischen den Punkten 0 und 1 als untere 

und obere Grenze. Da K, eine analytische Funktion ist, kann man den Weg zwischen Anfangs- 

und Endpunkt des Integrals beliebig wählen, ohne den Wert des Integrals dadurch zu verändern, 

solange man bei der Veränderung des Weges keine Singularität überschreitet, was vorausgesetzt 

werden soll. Das Element ds des Integrationsweges soll mit dem Linienelement ds, der Haupt- 

spannungslinie 1 den beliebig veränderlichen Winkel ß einschließen, so daß in der Gaußschen 
Ebene gilt: u 

dee dA TEE N ra: (19). 


Aus Gl. (14) folgt: 


1 1 
K(@) Fr PYı —D + im, — = Eike) de eat 

1 

en q ea (1 — 2]. —ia—Pp) .d 

-— aa +22% Ela, 24, e KJ a . (20). 
1 

= er = (a =]. iß 

= [len +202 jan 085 Ag 2 


Indem man auf beiden Seiten des Gleichheitszeichens Real- und Imaginärteil trennt, bekommt 
man 


1 1 
Ppı = R= (2% 0056 + 2 sin p)as—2 | 2: (& sinß + 2% c0sß)s un sil2daN, 
ö GES 085 a 05] 085 


Er 1 
00 [6 sinp + 22 eosp) as +2 | q- (cos — 52 sinp)as 28h. 
ö ö 


Wie schon erwähnt, kann der Integrationsweg zwischen dem Anfangspunkt 0 und dem Endpunkt1 
beliebig gewählt werden. Dieser Umstand dürfte den Hauptvorteil des neuen Auswerte-Ver- 
fahrens gegenüber den bisher bekannten darstellen. Man hat durch Wahl verschiedener Inte- 
grationswege, die alle zum gleichen Endpunkt 1 führen, die Möglichkeit einer sehr scharfen 
Kontrolle. Als Ausgangspunkt 0 wird man in Gl. (21a) wohl im allgemeinen einen Randpunkt 
wählen, wo p, bekannt ist. Um eine Integrationskonstante zu erhalten, ist es zweckmäßig, das 
Integral für einen Integrationsweg auszuwerten, für den nicht nur der Wert p, an der unteren 
Grenze, sondern auch p, an der oberen Grenze bekannt ist. 


$ 3. Ein Mittelwertsatz 
Wird in Gl. (21b) über einen geschlossenen Weg integriert, so verschwindet die linke Seite 
und es bleibt: 
$ (62 sinp + 510058) ds = 2-69: (2* cosß — 2° sinß)as 22 
05, 08, PER TR ea): 
Diese Gleichung gilt bis auf einen überall konstanten Wert von q, wie aus der Ableitung hervor- 
geht. Handelt es sich um eine geschlossene Hauptspannungslinie 1, über die integriert wird, so 
ist ds, statt ds und ß = 0 in Gl. (22) einzuführen, die sich damit vereinfacht zu: 


0% () 
264,1 ae (23). 


Diese Gleichung kann auch aus Gl. (15b), die man über eine geschlossene Hauptspannungslinie 1 
integriert, unmittelbar gewonnen werden, da 


ist. 


Gl. (24) kann als Mittelwertsatz gedeutet werden, wenn man den Mittelwert von E— AN 
für die geschlossene Hauptspannungslinie 1 folgendermaßen definiert: 2 
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Der Nenner 2 erklärt sich daraus, daß beim positiven Umlauf längs der Hauptspannungslinie 1 
der Winkel & sich um 2 x vermehrt: 


| | dan ER, Kann DEE ER N 
Der Mittelwertsatz läßt sich demnach so darstellen: 
1 .0q 


g= —$ —— ds 2 SMZ0N 
at, e 
Ist insbesondere die geschlossene Hauptspannungslinie 1 die Kontur eines lastfreien Loches, so 
ist an jeder Stelle der Lochkontur wegen o, = o, und 0, = 0 und Gl. (1) 


1 
a a a ee 


und damit läßt sich der Mittelwertsatz auch so ausdrücken: 


E u: 1 1 0 
29=-2p=-0=-, fan a=,_ = $1 
Lochkontur Lochkontur 


2.2. (29). 


Dieser Mittelwertsatz bezieht sich auf ein lastfreies Loch in beliebiger Größe in einem ebenen 
Spannungszustand. Er kann als Kontrolle eines spannungsoptischen Versuches dienen, bei dem 
das Isochromatenfeld an jeder Stelle des ebenen Spannungszustandes den Wert von q liefert, 
so daß längs der Kontur des Loches die unter dem Integral von Gl. (29) stehende Ableitung von g 
normal zum Lochrand als bekannt anzusehen ist. Besondere Vorsicht ist auf Grund unserer 
Erfahrungen, die wir bei der Kontrolle von Gl. (29) gemacht haben, in der Umgebung eines 
singulären Randpunktes, wo o, verschwindet, erforderlich, da hier 9g/os, unbestimmt wird. 


Bei dem hier abgeleiteten Mittelwertsatz sind die für den Lochrand definierten Mittelwerte 
von p und q gemäß Gl. (25) auf die Winkeländerungen da bezogen, die die Tangente an den 
Lochrand beim Fortschreiten längs des Lochrandes erfährt. Beim vollen Umlauf beträgt die 
Summe dieser Änderungen 2x. In einer vor Jahresfrist erschienenen Arbeit des Verfassers war 
ein Mittelwertsatz mitgeteilt worden, wobei der Mittelwert der Größen o,, p und q auf die Linien- 
elemente der Lochkontur bezogen wurden. Der an jener Stelle mit dieser Definition des Mittel- 
wertes ausgesprochene Mittelwertsatz ist nicht für jedes beliebig gestaltete Loch gültig, wie sich 
inzwischen herausgestellt hat 2). 


2) G. Sonntag, Integralsatz. Z. angew. Math. Mech. 34 (1934) S. 469—471. 
Eingegangen am 27. Dezember 1953. 


Die Einschließung von Eigenwerten hermitescher Matrizen 
beim Iterationsverfahren 


Von E. Kreyszig in Tübingen 


Im Anschluß an eine Theorievon Wielandt [1,2]') werden Einschließungsintervalle für die Eigen- 
werte einer gegebenen hermiteschen Matrix A auf Grund einer ebenfalls gegebenen Iterationsfolge von Vektoren 
%p = AR 5 = As, ermittelt. Abschnitt 1 bringt Ergebnisse von Wielandt, Abschnitt 2 
Beiträge zur Theorie des „‚Zweischrittverfahrens“ und Abschnitt 3 die numerische Durchführung der Einschlie- 
Bung beim Bin- und Zweischrittverfahren. 


Based on atheoryof Wielandt[1,2]!), inclusion intervals for the eigenvalues of a given Hermitian 
matrix A are found by aid of a likewise given iterative sequence of vectors x, m —=A % .- »%s = AS xy. Part 
gives theresultsof Wielandt, part 2 contains a contribution to the theory of the two-step-method, and part 3 
deals with the numerical performance of the process, üf the one-step- or the two-step-method are applied. 

En rapport avec une theorie de Wielandt [1,2]!) des intervalles limitants pour les valeurs propre 
d’une matrix Hermötienne donnee A sont trouves sur la base d’une sequence d’iterations donnee, de meme, de 
vecieurs &, % = A &y.. = As. Part 1donne des rösultats de Wielandt, part2 des contributions 
& la theorie de la methode de deux pas, et part 3 V’execution numerique de la limitation appliguee a la methode 
d’un ou de deux pas. 

Mcxons us Merona Bunannra [1, 2]1) ompenenamTca 3AKJTOUAINUE HHTEpBasst AH 
COÖCTBEHHEIX SHAYCHHH 3aNaHHOH 3PMHTOBOK MaTpulst A Ha OCHOBe TAKIKe MAHHON HTEPa- 


1) Laufende Nummer des Literaturverzeichnisses am Schluß der Arbeit. 
31* 


460 Kreyszig, Die Einschließung von Eigenwerten hermitescher Matrizen Bau 


IHOHHOH TOCHENOBATENBHOCTH BEKTOPOB 2, 44 = Ay, .-.; 1 — 4:0. B oTgene 1—M AaToH 

pesyssrarsı Bunangra, B orgere 2—M paspabatzıBaeTca TeopHuA ABYXINATOBOTO METOAA, 

& B 3—M UHCIIEHHO BBINONHAITCA MeTOABI OAHOTO H ABYX IIATOB. 

Gegeben sei eine hermitesche Matrix A und ein Ausgangsvektor 2, # 0. Aus der durch das 
Iterationsverfahren gelieferten Folge von Vektoren 


ie rer 


kann man Näherungen für die Eigenwerte },> 4, >...> 4, von A gewinnen. Wertvoll sind 
derartige Näherungen jedoch nur dann, wenn sich Schranken für die Abweichung von den exakten 
Werten angeben lassen. Dies geschieht mit Hilfe von Einschließungssätzen ?). Eine von W ie- 
landt [1, 2] entwickelte Einschließungstheorie gestattet es, die in die Iteration hinein- 
gesteckte Arbeit vollständig zu verwerten, indem alle Vektoren der Iterations- 
folge berücksichtigt werden können. Die vorliegende Untersuchung befaßt sich mit der prak- 
tischen Durchführung dieser Theorie 3). Als Grundlage werden in Abschnitt 1 die benötigten 
Definitionen und Sätze aus [1] und [2] zusammengestellt. Von der Auswertung zusätzlicher 
Vorkenntnisse, z. B. über eigenwertfreie Intervalle, etwa auf Grund bekannter positiver Definit- 
heit, wird abgesehen. Diese Fragen sollen in einer späteren Arbeit behandelt werden. 

Leitfaden für den praktischen Gebrauch: Einen allgemeinen Über- 
blick mit beschränkter Genauigkeit über die Lage und Länge der Einschließungsintervalle kann 
man sich nach Satz 4 mit Hilfe der nach Hilfssatz 6 zu konstruierenden Momentenellipse ver- 
schaffen. Satz 7’ liefert dieLänge 2d(y) des Einschließungsintervalles zu gegebenem Mittelpunkt y. 
Es gibt zwei Werte von y, y—=y, und y= y,, für die 2d(y) am kleinsten wird. y, und %s 
findet man durch Satz 8, praktisch brauchbare Näherungen für y, und y, nach Satz 10 und 
Schranken für d(y,) und d(y,) nach Satz 11. Wesentliche Vorteile hinsichtlich der Genauigkeit 
und des Arbeitsaufwandes der numerischen Rechnung bietet die in Abschnitt 2.3 behandelte 
Iteration mit Orthogonalisierung. 


1. Ergebnisse von Wielandt [l, 2] 


Definitionen: 1. Ein Spektrum von Eigenwerten A, >4, >--.->/, heißt 
verträglich mit den Vektoren %,, &ı. -. . , %,, wenn es eine hermitesche Matrix A mit den 
Eigenschaften 


n 
Ay, =, 2 gen. rn, air, , det (AN DEIN) 
i—1 
gibt ?). 

2. Eine Teilmenge F der reellen Zahlengeraden @ heißt eineEinschließungsmenge, 
wenn sie aus jedem mit 2, %; ,...,%, verträglichen Spektrum mindestens einen Eigenwert 
enthält. Insbesondere heißt F eine minimale Einschließungsmenge, wenn 
kein echter Teil von F Einschließungsmenge ist. — Eine Einschließungsmenge, die ein ab- 
geschlossenes Intervall auf der projektiven Geraden ist, heißt ein Einschließungs- 
interwa ll, 

3. Wir setzen 


= (2, u) = (4un Ar); 0, TE A ee 


und bezeichnen diese Skalarprodukte der Vektoren x, &,,...,2©,als Momente. 

m, ist bei hermiteschen Matrizen unabhängig von k. 

Die Berücksichtigung von 1,2,...,9% (<s—1), aufeinanderfolgenden Iterations- 
schritten, d.h. von 2,3,...,p9 + 1 aufeinanderfolgenden Vektoren &; der Iterationsfolge bei der 
Eigenwerteinschließung bezeichnen wirals Ein-, Zwei-,..., »Schrittverfahren. 


Satz1 [1]. Beim Einschrittverfahren gibt es zu jeder Zahl « genau zwei 
minimale Einschließungsintervalle, für die v Endpunkt ist; ihr anderer Endpunkt v;(w) ist durch 


vu) = (me i+ı U magi+2)|(mai U Mg i-+1) 5 = 0, 1, RE 2s—2 5 
gegeben. (Eines von diesen Intervallen enthält den uneigentlichen Punkt o). 
v0) = Oi, %i41) = Meirılmai, ae d,. 2 0 


®) Vgl.z.B.Collatz [3], S. 167, wobei jedoch nur ein Teil der in die Iteration gesteckten Arbeit ver- 
wertet wird. Vgl.auch Bückner[4]. 

) Herrn Prof. Dr. Wielandt danke ich herzlich für die Anregung zu dieser Arbeit und zahlreiche 
wertvolle Diskussionen und Hinweise, der Deutschen Forschungsgemeinschaft für die Unterstützung bei dieser 
Untersuchung. 


4) Statt ,„Eigenwert‘findet manin der Literatur auch die Bezeichnung ‚charakteristische Zahl“, 
z.B. bei Collatz [3]. 
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ist eine untere Schranke für den größten Eigenwert A, von A und zugleich die Mitte des absolut 
kürzesten minimalen Einschließungsintervalles; dessen Länge beträgt 2 ö(&;, &11), wobei 
6°(2;, +1) = me i+2|Mmg ı — 0% %+1) we 0, 1, ...), 28 RR 2 9 ist. 

v;(u) , oe und ö hängen von i, d.h. davon ab, welcher Iterationsschritt für die Ein- 
schließung herangezogen wird. 

Die Theorie von Wielandt gilt für beliebige endliche Schrittzahl p. Mit der Anzahl 
berücksichtigter Schritte wächst zwar die Güte der Einschließung, zugleich aber auch der 
Rechenaufwand. Es dürfte im allgemeinen nicht zweckmäßig sein, über das Zweischrittverfahren 
hinauszugehen. Dieses liefert, wie sich zeigen wird, wesentlich bessere Ergebnisse als das Ein- 
schrittverfahren, wobei sich die Rechenarbeit im Vergleich mit der Ausführung eines weiteren 
Iterationsschrittes bei Benutzung der im folgenden entwickelten Näherungsverfahren lohnt. — 
Wir formulieren deshalb die Ergebnisse von Wielandt nur für das Zweischrittverfahren. 

Ein aus einer Mindestanzahl von Eigenwerten bestehendes verträgliches Spektrum heißt 
ein Minimalspektrum. Beim Einschrittverfahren besteht dieses aus dem beliebig wähl- 
barem Wert u und dem zugehörigen Wert v;(u) nach Satz 1. 


Satz2. Beim Zweischrittverfahren werden alle minimalen Einschließungsinter- 
valle jeweils durch zwei auf der projektiven Geraden benachbarte Eigenwerte eines Minimal- 
spektrums .begrenzt. Dieses besteht aus den drei Eigenwerten A, A, und A,. Zu beliebig 
gegebenem A, erhält man A, und A, als Wurzeln von 


M MM Mm, 1 


_ |Mı Ma M; A, 
A m; m; m, Ai 
AED 

Geometrische Deutungen der Verträglichkeit stützen sich auf die 


4 0 et ya (). 


Definition 4. Ein Punkt mit den kartesischen Koordinaten 


n(a) = m,(a)/m,(a) und Z(a) = mz(a)/my(a) , oT To 70. BE FEZT, 
wobei 
m,(a) = (Ku Ar); re 0pl52, ee ala), 
die zu den Vektoren 
gehörenden Momente sind, heißt ein Momentenpunkt M(a). Die Menge aller M(a), 
die eine Ellipse bildet, heißt die Momentenellipse E. — In demselben Koordinaten- 
system heißt ein Punkt mit 7) = A, und &w)=% ein Eigenwertpunkt L(4,) 
Mit den Bezeichnungen von Satz 2 und Definition 4 gilt 
Satz 3. Die Strecke L(A,) L(A,), i# k, also die Verbindung der Eigenwertpunkte zweier zu 
demselben Minimalspektrum gehörenden Eigenwerte, berührt die Momentenellipse im Momenten- 
punkt M(—- Ar’), k#lz#i. 
Satz 4. Für jedes verträgliche Spektrum 4, >%>:---> A„liegt dieMomentenellipse notwendig 
in dem konvexen n-Eck mit den Ecken L(A,). 
Nach dieser Zusammenstellung der verwendeten Hilfsmittel wenden wir uns der Behandlung 
des eingangs formulierten Problems zu. 


2. Beiträge zur Theorie des Zweischrittverfahrens 
2.1 Geometrische Methoden 

Die geometrischen Methoden knüpfen an die Definition 4 an. Sie vermitteln einen all- 
gemeinen Überblick über die Einschließungsintervalle, sind also gut geeignet, sich mit dem Zwei- 
schrittverfahren vertraut zu machen. Ihre Genauigkeit ist jedoch beschränkt, deshalb muß eine 
präzise Bestimmung der Einschließungsintervalle analytisch erfolgen, wie in den folgenden Ab- 
schnitten gezeigt wird. 

Zuerst geben wir ein Verfahren zur Konstruktion der Momentenellipse an. 


Definition 5. Die Unterdeterminanten B;, zu den Elementen b;; der Determinante 
N, My. Me 
B=|b.l= |m, m, m; 
Ma My; My 
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werden durch 


ke kz ky 
(Bir m k; k; kz 
k ka ka 
bezeichnet, also ausführlich 
| ma My My » 1m Img TeRS Sup _ | m, m; 
Im, ma|’ M, Ma” Mg Mg Mg, My M; My 
und 
ko = Mg Ms 
M; My a 
Sie heißen Momentenkombinationen. — Ihr Index ist gleich der Summe der 


Indizes der Produkte m, m, gewählt. 
Hilfssatz 5. Die Momentenellipse Z ist durch 
M Ms; 1 |mı m, 2 
mı Mm n |m m, 2n 
Mm Mm € Im mı 27 ns kö—kıntk; 2(E—knt ke) el) 
My My 9): My Ma 1 2(k,&— kn + kı) k&—kıntk; 
mM Mm 2n m, m n 
m; m; 2 m m, & 
gegeben. 
Beweis. Aus (3) folgt 
m,(a) = m, + 2am,+ı + a? mı 12, vr 0. 1.020. a (3), 
und hieraus wegen (2) für einen Punkt (n, &) der Momentenellipse E 
am; n — ms) +2 um n — m) mn —mı —0, 
am, & — mı) + 2 am, &E— m) Hm &E— m, — 0 


Bekanntlich haben nun zwei Gleichungen a2? +2ßzx +y=0 und A®+2Bz+(0=0 
eine gemeinsame Wurzel, wenn ihre Resultante 


Acl® Aa| |Bß 
R= 0, * BB le, 
verschwindet. Im vorliegenden Fall wird 
Mg Ma ıP Immo m, 11 |mı m, 1 
R= m, m; n) —4 m, m, n) |m, m; n) = 0 
m; m & m; m; E| |m; m, E 


Dies ist mit der Behauptung identisch. 


Die praktische Konstruktion der Momentenellipse E erfolgt am einfachsten punktweise, 
gestützt auf den 


Hilfssatz 6. Die Abszissen n(@) der Momentenpunkte M(a) liegen zwischen den Wurzeln R, 
und R,< R, der Gleichung 

MA)=k,A—k,A+k=0 (vgl. Def.4.und5) ... 2... (5). 

Zu den Abszissen n gehören die Parameterwerte 
an) = — 95" (9 + 91 9) 
mit 
Gi—=Mm —NMi-ı; i—=1,2,3; 

und zu diesen die Ordinaten &£ nach Definition 4. 
Beweis. Aus Satz 3 folgt, daß R, und R, zusammen mit A= o ein Minimalspektrum 


bilden. Also liegen nach Satz 4 alle Abszissenwerte 7 von E in dem abgeschlossenen Intervall 
[R,, R.]. Die Beziehung für a(n) folgt unmittelbar aus der Gl. (4). 


ee a 


Da ngew; Math. Mech. Kreyszig, Die Einschließung von Eigenwerten hermitescher Matrizen 463 


Eine rationelle Näherungsmethode zur Bestimmung von E besteht 
darin, daß man außer R, und R, noch ein oder zwei geeignete Minimalspektren nach Satz 2 
berechnet. Wegen Satz 3 hat man dann ein Tangentensechs- oder -neuneck, in das man E bequem 
eintragen kann. 

' Die geometrische Darstellung der Einschließungsintervalle kann nach Wielandt noch 
übersichtlicher gestaltet werden, indem man das Gebiet £{> n? nach Definition 4 vermöge 


yn g— v.—n? GE © Mo. 10 «Br een en, (6) 


in das Poincar&sche Modell der hyperbolischen Ebene transformiert, wobei die Parabel 
&=n? mit den daraufliegenden Eigenwertpunkten in die y-Achse und alle Parabelsehnen in 
Halbkreise mit dem Mittelpunkt auf der y-Achse übergehen. Dann gilt 


Satz 7. E’ sei das Bild der Momentenellipse E bei der Transformation (6). Dann erhält man 
die Enden eines minimalen Einschließungsintervalles zu gegebenem Intervallmittelpunkt y als 
Schnittpunkte der y-Achse mit demjenigen Kreis um den Punkt y, der E’ berührt. — Ferner 
durchläuft der Intervallmittelpunkt y zweimal monoton die Abszisse, wenn der Berührpunkt 
des zu y gehörenden Berührkreises von E’ einmal monoton die Kurve E’ durchläuft. 


Beweis. Der innerhalb der Parabel & = n? gelegene Teil der Tangenten an E geht durch die 
Transformation (6) in Berührhalbkreise von #’ mit dem Mittelpunkt auf der y-Achse über. 
Aus Satz 4 folgt damit der erste Teil der Behauptung. Die monotone Abhängigkeit der Berühr- 
kreismittelpunkte von den Berührungspunkten ergibt sich durch einen Vergleich der Tangenten 
an E mit ihren Bildern bei der Transformation (6). 

So kann man in der yz-Ebene die Lage und Länge der Einschließungsintervalle unmittelbar 
übersehen, wie dies Bild 2 zeigt. Satz 7 umfaßt als Sonderfall die geometrische Konstruktion 
der Weinsteinschen Einschließungsintervalle des Einschrittverfahrens ([1], Seite 240). 


2.2 Analytische Methoden: Minimale Einschließungsintervalle zu gegebenen Mittelpunkten, 
kürzeste minimale Einschließungsintervalle 
Die analytische Formulierung des Satzes 7 folgt aus Satz 3; es gilt 


Satz 7‘. Die halbe Länge d(y) eines minimalen Einschließungsintervalles zu gegebenem 
Mittelpunkt y ist die Quadratwurzel aus der kleineren Nullstelle der quadratischen Gleichung 


Mo mı Mg 1 

Mm; My Ms y+d | k,d“ - 
m m m YHd| Zap y ti) 2% 
1 y-dAg-D 0 | up ah Hit Lkyt ht. 


Ferner liegen im Intervall [y — d,(y) , y + d.(y)], wobei d3(y) die andere Nullstelle von (7) 
ist, zwei Eigenwerte jedes verträglichen Spektrums, weil nach Satz 3 jedes Minimalspektrum 
beim Zweischrittverfahren drei Einschließungsintervalle liefert. 

Von besonderem Interesse sind die kürzesten minimalen Einschließungsintervalle d.h. die- 
jenigen, deren Länge 2 d(y), absolut genommen, einen kleinstmöglichen Wert hat. Deren Mitten 
gibt 
Satz 8. Die größte Wurzel, y,, und die kleinste Wurzel, y,, des Polynoms 

4 
2 tl a Der alu Sue SE ee (8) 


n—0 


a=(k—Ak,k)k,, 
g=—hztkek(6ky—k,)+2kik,, 

g=k(kk —k)— 2Akzkzkz + k,(k/2 — k,)? — 4 k, ki, 
sk (ik —krlA+kyke) + 2kykık,, 

= (1/4) (k;k; (k, — 2ky— köke — k,k;) ’ 


sind Mittelpunkte von kürzesten minimalen Einschließungsintervallen. Die anderen beiden 
Wurzeln, y, und %,, sind Mitten längster minimaler Einschließungsintervalle. 


mit 


Beweis. (8) folgt aus (7), indem man die Ableitung von d nach y null setzt. Die übrigen 
Behauptungen werden zusammen mit Satz 10 bewiesen. 

Das zu y, gehörige Intervall F(y,) wollen wir als „zweitkürzestes“ bezeichnen, es 
ist kürzer als Intervalle mit Mittelpunkten y in einer Umgebung von 4,; trivialerweise existieren 
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beliebig viele Intervalle, die kürzer als F(y,) sind, mit Mittelpunkten in ein geeigneten Um- 
ebung von %,- 

: Zu einfacheren Polynom, als es (8) ist, gelangt man, indem man den pralstäsch 

wichtigen Intervallmittelpunkt y durch eine theoretisch angemessenere Variable, A, ersetzt: 


i i ein Mini t das minimale 
Hilfssatz 9. Bilden A, A, und A, ein Minimalspektrum nach Satz 2, so hat da 
Einschließungsintervall 1A, A, extremale Länge, wenn für A, der Wert einer der vier Wurzeln 
des Polynoms . 


Mm Mı Ma mM Mı Ma 
p(A)= Im, m; m (my ?— 2m, At m) —2|mı m, mM; (m #2 — 2m, iA+m)=0 (0 
11 22. 2998 1v 9b 
gewählt wird. > 
Beweis. (9) folgt aus (1) durch geeignete Differentiation unter Vereinfachung der Deter- 
minanten aus den k,?—=2,...,6 (vgl. Def. 5) mit Hilfe des Sylvester schen Determinanten- 
satzes. Setzt man die Wurzeln von (9) in (1) ein, so erhält man die Einschließungsintervalle 


extremer Länge in Übereinstimmung mit (7). : \ ; 
(8) hat gegenüber (9) den Vorteil, daß man bequem gute Näherungsaussagen über die 
Wurzeln %,, Y5, y, und y, machen kann. Hierzu dient 


Satz 10. Die Mitten kürzester minimaler Einschließungsintervalle nach Satz 8 sind näherungs- 


weise durch 
Yı.=,R, . und Iys = har MV@L-Filtssatz 077, 2 m (10) 


und in besserer Näherung durch 
Mer ed ya Eee EEE (11) 
mit 
e* = e—= D"?{(3 k2]4 k2 — kylk,) D? + c, + c5k,]2 k,— 3/4 DEN , D=R,—R,, 
gegeben, wobei c, und c, durch Gl. (8) definiert ist. Ferner gilt (vgl. Satz 8) 


rk ek een a De (12) 
und verbessert 


mit 
y=kslk; + 03/2 DE. 
Ist ferner y, #£ %s, so ist 4 
&Se und d(y)Zd(y,) für yZ0. 


Vorbemerkung: eundy sind Differenzen großer Zahlen, ihre prinzipiell einfache Berech- 
nung muß deshalb hochstellig, etwa mit 6 bis 8 Stellen durchgeführt werden. 


Beweis des Satzes 10 und Schluß des Beweises des Satzes. 


E’ hat nach Satz 7 für y= Y,, Ya Y5 Y, extremale Ordinatenwerte z(y). Die Transfor- 
mation (6) führt die große Halbachse von Z in ein Kreisbogenstück über, dessen Mittelpunkt 


rt 
i. a. von der Abszisse des Mittelpunktes von E, 
Yy= ks]2k,, 


die sich aus Hilfssatz 6 ergibt, verschieden ist. Kleinste Werte z(y) liegen desto näher unterhalb 
von R, bzw. oberhalb von R, (vgl. Hilfssatz 6), je gestreckter Z ist, d. h., wie aus dem Entwick- 
lungssatz folgt, je mehr die Beteiligung zweier Eigenwerte an % gegenüber den anderen 
überwiegt. Dabei heißt ein Eigenwert 4; beteiligt an x, wenn das Skalarprodukt der zu 
A; gehörenden Eigenlösung mit x, nicht verschwindet. Zwei Wurzeln des Polynoms (8) — wir 
bezeichnen sie mit y, und y, —- sind also Mitten kürzester Einschließungsintervalle; so 
ergibt sich die Berechtigung für die Näherungsausdrücke (10) und (11);es ist e> O und e® > 0. — 
Ein längstes Intervall muß das Unendliche enthalten; seine Mitte bezeichnen wir mit Y. Da E' 
eine für alle Werte y#£ R, und y £ R, differenzierbare Jordankurve ist, muß sich zwischen 
y, und y, ein Maximum des der Abszisse zugekehrten Teiles der Kurve E’ befinden; die zu diesem 
Maximum gehörige Abszisse %, ist wegen Satz 3 Mitte eines (relativ zu Intervallen mit benach- 
barten Mitten) längsten Einschließungsintervalles. Aus dem Vergleich von E mit E’ folgt, daß 
Y, und y, in einer Umgebung von y —= m liegen, die sich im Falle der Entartung von E in eine 
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Strecke auf den Punkt y = m zusammenzieht. Hieraus folgen die Gleichungen (12) und (13). Die 
Werte e und y liefert der Vietasche Wurzelsatz, indem &* — e gesetzt wird, unter Vernach- 
lässigung von Gliedern mit e?. Für y = 0 ist „= m, dann ist E’ axialsymmetrisch bezüglich 
der Geraden y—= m und damit gilt wegen Satz 7 die Beziehung d(y,) = d(y,) sowie &* —e. 
Ist y>0, also m < y, so folgt aus den Eigenschaften der Transformation (6), daß 
E(y +w <Eky—w, (ytw<R, y—w>R,), und e* > e ist, daraus ergibt sich 
d(y,) > d(y,). Für y < 0 ist es umgekehrt. D. h., ob das absolut kürzeste minimale Einschlie- 
Bungsintervall zu 4, oder zu %, gehört, erkennt man aus dem Vorzeichen von y. 

Beim Einschrittverfahren bildet og mit A= oo ein Minimalspektrum und ist Mitte des 
kürzesten Einschließungsintervalles (vgl. Satz 1). Beim Zweischrittverfahren bilden R, und R, 
mit A= o ein Minimalspektrum und sind asymptotisch, d.h. bei Entartung der Momenten- 
ellipse E in eine Strecke, Mitten kürzester Einschließungsintervalle, wie soeben gezeigt wurde. 

Nicht nur für die Lagen kürzester Intervalle, sondern auch für deren Längen kann man 
brauchbare Näherungen angeben: 


Satz 11. Für die halben Längen d(y,) und d(y,) der kürzesten minimalen Einschließungs- 
intervalle existieren Schranken, die nur von den Näherungswerten für y, und y, nach Satz 10 
abhängen; es ist 


U(R;) > d(y) > d(R;) (1 — 4 &/K RP)? = d(R,) — 2e2/d(R)) , nl, 2a 


Beweis. Die linke Ungleichung ist zufolge Satz 7 und 10 trivial. Es ist |R; — y| S2.:=- 12 
und der Faktor2 kann durch keinen kleineren ersetzt werden. Ein Kreis durch den Punkt P 
mit den Koordinaten y— y, und z—= d(y,) mit dem Mittelpunkt (y= R,, 2 = 0) schneidet die 
Kurve E’, weil es zu P wie zu jedem Punkt von Z’ genau einen Berührungskreis gibt und dieser 
für P den Mittelpunkt (y—= y,2=0) hat. Also gilt d(y,)? + 2Qe)?>d(R,)” und das Analoge 
für y, und R,. Hieraus folgt die rechte Ungleichung. 

Satz 11 liefert desto bessere Schranken, je stärker &,, x, und x, näherungsweise proportional 
zueinander sind. 

Hat man Lage und Länge des kürzesten minimalen Einschließungsintervalles, z. B. nach 
Satz 10 und 7’, festgestellt, so kann man die Länge des zweitkürzesten minimalen Einschließungs- 
intervalles (vgl. Satz 8) roh abschätzen vermöge 


Hilfssatz 12. Ist die Mitte y, des kürzesten minimalen Einschließungsintervalles von der 
Mitte y, des zweitkürzesten (vgl. Satz 8) verschieden, dann gilt 


a > EN NE UEREREER (15) 


und zwar desto genauer, je geringer die Differenz A der beiden zun = k;/2 k, gehörenden Werte 
sind. R, R, und £ sind durch Hilfssatz 6, k, und k, durch Definition 5 und y durch (13) gegeben. 


Beweis. R,—R, ist nach Hilfssatz 6 und Satz 7 die Abszissenerstreckung von E bzw. E’. Die 
Transformation (6) führt die große Achse von E in ein Kreisbogenstück über; den zugehörigen 
Kreis bezeichnen wir mit K. z, bzw. 2, seien die Ordinaten der Schnittpunkte von X mit den zur 
z-Achse parallelen Tangenten an E’, die durch y= R, und y= R, gegeben sind. Elementar- 
geometrisch findet man 2 = 27 +2y(R, — R,). n= ks]2 k, ist die Abszisse des Mittelpunktes 
von E, also ist A ein Maß für die ‚Dicke‘ von E. Entartet E in eine Strecke, also E’ in ein Kreis- 
bogenstück, d.h. ist A = 0, so wird 2,— d(R,) und 2, = d(R,), wie behauptet. 


2.3 Iteration mit Orthogonalisierung 
Der praktische Wert der soeben entwickelten Methoden hängt zu einem guten Teil von der 

Genauigkeit ab, mit der man die durch Definition 5 gegebenen Momentenkombinationen k,, 
Is ‚kg, ermittelt. Diese k; sind Differenzen großer Zahlen. Wie aus dem Entwicklungs- 
satz (siehe z.B. [3]) folgt, werden die k; klein bzw. verschwinden, wenn z,, #, und x, näherungs- 
weise bzw. genau proportional zueinander sind. Deshalb ist es zweckmäßig, das Iterationsver- 
fahren so anzulegen, daß diese Nachteile wegfallen. Zum Ziele führt eine Orthogonalisierung 
nach jedem Iterationsschritt und bezüglich aller vorhergehenden Vektoren der Folge. Dabei 
können alle in dem Iterationsprozeß infolge der Orthogonalisierung zusätzlich auftretenden 
Rechengrößen bei der Bestimmung der k; wiederverwendet werden. Wie man durch Rechnung 
zeigt, gilt 
Satz 13. Die Vektoren u 

2, = A — 410% 


und 
2, — A2ı — Ay1?ı — 490% 


> 
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sind für | | 
a = mm, = (2 A) /@o 20) » 
a9 — Helm = (2ı21)/ 20%) 
Ag = Holz = (2 Azı)/2ı 2) 
orthogonal zueinander. Es gilt 

a) für die Momentenkombinationen nach Definition 5 


k,= Mole» k;— 09 + Aa» ka = Ayo Agı — Op » ’ 
k=b+ (ao + 9), k;—= 016 + (Gio@zı — 40) (dio + Azı) BEER TE) 
ka — b (ao + 0) + Asodzı (aioAzı — 2) + A 
mit 
b=uli, und. w=l22) 
b) für die Momente nach Definition 3 
My) in; Mr AH) = u: 
M—U;+Qpokı, M—Myt 20 He + Golı> RER, Di (17), 
mM + (20 + Azı) Us + Bao + Aa) He + ok 
und 


c) für die durch Satz 1 und Hilfssatz 6 eingeführten Größen 
OR 21) = Ayo » °(R %ı) = An 


(die Argumentangabe kennzeichnet den Iterationsschritt) und 


R,2 = (1/2) (ao + Qsı + Yan — 0, +4 PN) . 


Man geht also von &, = 2, aus, gewinnt durch den ersten Iterationsschritt Az,, orthogonali- 
siert dann, d.h., bildet z,, iteriert und orthogonalisiert den dabei erhaltenen Vektor Az, noch- 
mals, d.h. bildet 2, ebenfalls nach Satz 13. Die Formeln (17) zeigen, daß es sich nicht lohnt, 
zur Bestimmung der m; etwa zusätzlich den Iterationsschritt &, — x,—= A x, auszuführen. — 
Einen Teil der Größen von Satz 13 benutzt aus Karush ([5], Satz 1). 


3. Numerische Durchführung der Einschließung beim Ein- und Zweischrittverfahren 

Die numerische Durchführung der vorstehend entwickelten Methoden wird nun für Matrizen, 
gewöhnliche Differentialgleichungen und Integralgleichungen an je einem Beispiel gezeigt. Durch 
den Vergleich der numerischen Ergebnisse erkennt man die Überlegenheit des Zweischrittver- 
fahrens über das Einschrittverfahren. Die Reihenfolge der Überlegungen entspricht der Auf- 
einanderfolge der Sätze im vorhergehenden Kapitel 2. 


3.1 Numerische Durchführung der Einschließungbeieiner Matrix 
Für die Eigenwerte einer gegebenen symmetrischen Matrix 
5,68 3,24 —1,4 —0,6 —0,2 

3,24 7,57 1,05 0,45 0,15 

A=|—1A 1,05 3.25..,.0,25 0,75 

—0,6 0,45 025.032 1,75 

—0,2 0,15 0,78. 11475 3,25 
seien Einschließungsaussagen zu machen. (Diese Matrix wurde zu einem Spektrum konstruiert, 
wie es in praktischen Fällen etwa aufzutreten pflegt; es ist nämlich exakt /, = 10, A, — 6, 
A;—4, ,—=2 und 4,—=1.) Liegen irgendwelche Anhaltspunkte für die Lage der Eigenwerte 
nicht vor (und von derartigen „speziellen Vorkenntnissen‘ soll, wie eingangs erwähnt, hier kein 


Gebrauch gemacht werden, sondern erst in einer späteren Untersuchung) so können wir als Aus- 
gangsvektor für die Iteration z. B. 


%, = (7, 12, 4, 5, 6) 


mit den ganzzahlig gemachten Zeilensummen von A als Komponenten wählen und erhalten 
eine Iterationsfolge x, %,, 2, an der, wie sich zeigen wird, die Eigenlösung e, zu A, beteiligt 
ist, jedochnoch nicht zu sehr überwiegt, so daß die Allgemeinheit gewahrt bleibt und dieKurven Z 


und E’nicht übermäßig gestreckt sind, wie manausden maßstäblichen Zeichnungen (Bild1 und2) 
ersieht. 
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A. Einschrittverfahren: Nach Satz1 ergeben sich für A, —= 10 die unteren 
Schranken o(&, &,) = 8,715 beim ersten Schritt und o(z,, %,) = 9,505 beim zweiten. Diese 
Werte sind zugleich die Mitten kürzester minimaler Einschließungsintervalle von einer Länge 


[l 


Bild1. Beispieli, Momentenellipse Z im 
konvexen Fünfeck der Eigenwertpunkte 


Bild 2. Beispiell, Transformierte Kurve E ? 
im Halbkreisbogenfünfeck und kürzeste mini- 
male Einschließungsintervalle beim Ein- 


schrittverfahren (1. Schritt: Mittelpunkt 1 
0%, %,) 2. Schritt: Mittelpunkt o(z,, 2;)) 
und beim Zweischrittverfahren (Mittel- 
punkt y) een | 


2 ö(& &ı) = 2° 2,058 bzw. beim zweiten Schritt 2 ö(z,, %5) = 2 1,348, so daß man als kür- 
zestes Intervall 
beim ersten Schritt 


6,657 <A<< 10,774 
und beim zweiten Schritt 
8,157 <A<< 10,853 
findet. 
B. Zweischrittverfahren: Die Mitten kürzester Einschließungsintervalle sind 
die Wurzeln y, und y, der Gl. (8), für die nach Satz 10 die Näherungen 
Yyı = .R, = 9,878 und Ya Bs==5,.070 
und die besseren Näherungen 
Yyı = R, — e = 9,862 und Ya = R,; + e* = 5,086 
gelten. Ferner findet man mittels Satz 10 für die Mitten y, und y, längster minimaler Intervalle 
Ya = Ya = k;]2 k, = 7,47 und verbessert y; = %, = ka[2 k, — y = 7,306. Die mit der angegebenen 
Anzahl von Ziffern genauen Werte der Wurzeln von (8), wie sie z.B. das Newtonsche Ver- 
fahren liefert, sind 
Yı = 9,8634, Y; = 9,089 und = 421,305. 


Es lohnt sich also, e und y zu berücksichtigen, wie die folgende Tabelle zeigt: 


Abweichung vom exakten Wert | Yı | Ya |% und %, 
Rohe Näherung 1,50/o0 40/00 2,50), 
Verbesserte Näherung 0,1%/g0 | 0,6%/g0 0,5%/ 0 
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Für die Längen 2d(y,) bzw. 2d(y,) der kürzesten minimalen Einschließungsintervalle liefert 
Satz 11 die Schranken 


d(R,) = 0,877 12 d(y,) > 0,876 6 = d(R,) — 2 82/d(R)) 


d(R,) = 1,513 42 d(y,) > 1,513 1 = d(R,) — 2 e?/d(R,) 
und Satz 7 die in der angegebenen Stellenzahl exakten Werte 
d(y,) = 0,8771 bzw. d(y,)= 1,5134, 


so daß 8,9863 <A < 10,7405 das kürzeste und 3,5758 <A <s 6,6026 das „zweitkürzeste‘‘ 
(vgl. Satz 8) minimale Einschließungsintervall ist. Diese beiden Intervalle sind punktfremd 
zueinander, schließen also zwei verschiedene Eigenwerte von A ein. Die Näherungen für die Lage 
und Länge dieser beiden Intervalle nach Satz 10 und 11 fallen ungenauer als bei diesem Beispiel 
aus, wenn die Vektoren der Iterationsfolge in geringerem Maße annähernd proportional zuein- 
ander sind; im anderen Falle werden die Näherungen noch besser als hier. — R, nähert A, bis 
auf 1,2°/, an, R, den dem Betrage nach zweitgrößten an x, beteiligten Eigenwert A, = 6 bis auf 
15°%/,. — Aus dem nach Hilfssatz 6 und Satz 7 gezeichneten Original des Bildes 2 (Größe etwa 
20x30 cm) ergeben sich die Werte y, = 9,8, %; = 5,0, d(y,) = 0,9 und d(y,) = 1,6 mit einer 
Abweichung von der Größenordnung 1°/, von den errechneten Werten. 


und 


3.2 Numerische Durchführung der Einschließungbeieiner 
gewöhnlichen Differentialgleichung 


Im vorangehenden Abschnitt haben wir alle wesentlichen Fragen ausführlich behandelt und 
fassen uns deshalb hier kürzer. Die Differentialgleichung 


"+iy=0, y0=0, WU) ryd—0 


wurde von Gollatz ([3], Seite 163), ausgehend von der Funktion z,= 1, iterativ behandelt. 
Wir werten die ersten beiden der angegebenen drei Schritte aus. Unserem m; entspricht die 
Größe a;, unserem A die Größe 1/A bei Collatz. 


A. Einschrittverfahren: Die kürzesten minimalen Einschließungsintervalle 


sind 
| (Mitte) | (halbe Länge) | (Intervallenden) 
beim 1. Schritt 0(% %ı) = 0,208 ö(z. 2), 0,081 0,127, 0,289 
beim 2. Schritt 0(X1, %) = 0,2426 Ö(&1, %) = 0,0091 0,2335, - 0,2517 


B. Zweischrittverfahren: Die Mitten kürzester minimaler Einschließungs- 
intervalle sind nach Satz 10 näherungsweise 


Yı = R, = 0,242 750 3 und % = R, = 0,016 648 2 
und verbessert 
Yyı = Rı— e= 0,2427493 und %=R,-+ e* = 0,016 649 18. 
Die in der angegebenen Ziffernzahl genauen Werte sind 
Yy, = 0,242 749 1 und Y. = 0,016 649 19. 
Die relativen Abweichungen der Näherungen für y, und %, betragen 
(Rı — yıllyı = 410%, (IE y)lyı= 10% 
(ya — R;)/ya = 6.105, (ya— Rz — e*)/y, = 6.1077. 
era der angegebenen Ziffernanzahl genauen Werte der zugehörigen Intervallängen nach 
d(y,) = 0,006 719345 und d(y,) = 0,015 888 893 , 
ergeben sich auch bereits durch die Schranken nach Satz 11, da 
U(R,) — d(y,)<28/dR,) = 3-10-1% und UR,) — d(y;) <2 &2/d(R,) = 1, 2- 10-10 


er — R, nähert A, = 0,242 96 bis auf rund 0,1°/,an; R, ist als Näherung für A, völlig unbrauch- 
ar. 


= 
. 
2 


us. m 


\ 
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33 Numerische Durchführung der Einschließung bei einer 
Integralgleichung 


Schließlich behandeln wir die Integralgleichung einseitig eingespannter transversal schwin- 
gender Stäbe konstanten Querschnittes 


1 | wi (ER — 7/6 (>) 
SEEN z)dr—Ac) mit KE, = rn eos (E<9), 


ausgehend von u —=1. 


A. Einschrittverfahren: Die kürzesten Intervalle beim 1. bzw. 2. Schritt, 
0,0112 <A < 0,0888 bzw. 0,0798 <A < 0,0820 haben die Mitte o(&, x.) = 0,05 bzw. 
0(%1, &5) — 0,080 875, diese sind untere Schranken für A, = 0,080 89068, das sich aus 
w, = 1,875 10407 + 2-10-8 als größter Wurzel 5) der Gleichung cos w &o| w= — 1 ergibt, 
es ist A, = wit. 


B. Zweischrittverfahren: Für die Mitten 
Yı = 0,08088305 bzw. %s= 0,001 029 929 

der kürzesten minimalen Einschließungsintervalle 

0,080 103 <A < 0,081 663 bzw. 0,000 048 <A < 0,002 012 
nach Satz 7’ und 8 ergibt Satz 10 die Näherungen 

Yı = R, = 0,080 88306 bzw. %, = R, = 0,001 029 926. 

Deren relative Abweichungen 

(RR, — y)IRı=10” bzw. (y»— R)/R,=3.10* 
werden durch die Berücksichtigung von e = 4. 10”? auf 6. 10? bzw. auf 10° gesenkt. 


Der Vergleich der numerischen Werte des Einschrittverfahrens mit denen des Zweischritt- 
verfahrens zeigt die Überlegenheit des letzteren hinsichtlich der Schärfe der Einschließungen 
und des Umfanges möglicher Aussagen. Ein entsprechendes Bild ergibt sich, wenn man zusätz- 
liche Vorkenntnisse, z.B. über positive Definitheit berücksichtigt, wie dies in einer späteren 
Untersuchung geschehen soll. 
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KLEINE MITTEILUNGEN 


Ein Integralsatz der ebenen Elastizitätstheorie. ses verschwindet die Spannungssumme (co, + oy) des 
Eingangs möchte ich nicht versäumen meinem hoch- durch die eingekreiste Belastung erzeugten Spannungs- 

verehrten Lehrer Herrn Prof. Dr. L. Föpp| für die zustandes. Siehe Abb. 1; formelmäßig ausgedrückt: 

Anregung zu diesem Aufsatz zu danken. a $ EEE (). 
Es wird ein Integralsatz bewiesen und mit Hilfe des- 

selben der Föpplsche Mittelwertsatz!) für Kreis- 

löcher bestätigt. Weiter wird auf ein Versehen in der 

Beweisführung dieses Mittelwertsatzes aufmerksam ge- 

macht und an Hand eines Beispieles gezeigt, daß er 

nicht allgemein gültig ist. 
Der Integralsatz, der nach Kenntnis des Verfassers 

noch nicht ausgesprochen wurde, lautet: 
Eine unendliche Scheibe von überall gleicher Stärke 

sei in ihrer Ebene durch beliebige endliche Kräfte und 

Momente belastet; über den Weg eines beliebigen Krei- 


1) L.Föppl: „Ein Mittelwertsatz derebenen Elastizitätstheorie“. 
Aus den Sitzungsberichten der Bayrischen Akademie der Wissen- 
schaften, Mathematisch-naturwissenschaftliche Klasse 1952 Nr. 17, 
S. 215 bis 217. Bild 1 
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Beweis: Man denke sich die Scheibe längs des 
kreisförmigen Integrationsweges ausgeschnitten und 
belaste den Lochrand mit o, und z,; entsprechend den 
dort vorhanden gewesenen Spannungen. Für die äußere 
Scheibe einschließlich des Lochrandes hat sich damit 
nichts geändert. Für eine konstante Normalbelastung 
des Lochrandes ist überall o, +0; = 0 und für eine 
konstante Schubbelastung ist o,—= 0: = 0°). Daraus 
folgt, daß auch für eine beliebige Belastung normal 
oder tangential zum kreisförmigen Lochrand das Inte- 
gral Gl. (1) verschwinden muß, denn wäre dieses Inte- 
gral für eine Teilbelastung über a da des Lochrandes 


TR 


KG 


Bild 2 


nicht Null, dann könnte auch die Summierung dieser 
Integrale für eine gleichmäßig über den Lochrand ver- 
teilte Belastung nicht verschwinden. Mit anderen Wor- 
ten: DasIntegral Gl. (1)ist unabhängig von der An- 
griffsstelle einer Teillast p-a da an einem kreisförmi- 
gen Lochrand. Wenn mehrere derartige Teillasten mit 
gleicher Kraft und im gleichen Sinne (z. B. alle radial) 
am Lochrand angreifen, dann ist der obige von einer 
Teillast herrührende Integralwert mit der Anzahl 
dieser gleichartigen Belastungsteile zu multipli- 
zieren. Wird zu einer gleichmäßigen Lastverteilung 
übergegangen, dann verschwindet obiges Integral 
Gl. (1) (weil überallo, = —o;). Wenn aber eine Summe 
von Integralwerten gleichen Vorzeichens verschwin- 
det, dann ist das Integral jedes Teiles dieser Summe, 
d.h. für jede beliebige Teillast am Lochrand gleich 
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Diese Beweisführung läßt erkennen, daß obiger 
Integralsatz unverändert gilt, wenn der Integrations- 
weg kreisförmig ein beliebiges Loch umschließt, 
s. Abb. 2. r 

Es werde jetzt untersucht, wie sich hieraus ein Mittel- 
wertsatz ableiten läßt. Betrachten wir eine lochfreie 
Scheibe, die an ihren Rändern und auch in ihrem 
Inneren in ihrer Ebene belastet ist (Abb. 3;) dafür sei 
unter Voraussetzung elastischer Beanspruchung die 
harmonische Ortsfunktion der Spannungssumme mit 


q(2, y) = (02 + 0y)o 


bezeichnet. Der arithmetische Mittelwert der Span- 
nungssumme über einen Kreis auf dessen Rand und in 
dessen Inneren keine Singularitäten liegen, ist: 


ei” 248). 
= (8) 
Nach dem Mittelwertsatz der Potentialtheorie *) ist 
dieser Mittelwert 7 gleich dem Funktionswert der 
Spannungssumme im Mittelpunkt M dieses Kreises. 

Wird jetzt innerhalb des kreisförmigen Integrations- 
weges ein Loch in die Scheibe geschnitten, ohne die 
äußeren Kräfte zu ändern, dann ist der ursprünglichen 
Funktion g(x,y) eine Störfunktion p(x, y) zu über- 
lagern. 

Die resultierende Spannungssumme sei mit dem 
Index 1 bezeichnet, dann folgt die Gleichung 


(0. + oy)ı= (04 0y)o+ Plz, y)= 4la, y)+ plz, y) (@). 


Unter der Störfunktion o(z,y) ist die Spannungs- 
summe (0, + 0y)= p(x, y) verstanden, die aus der 
Belastung des Lochrandes herrührt, die man nach dem 
Ausschneiden des Loches aus der belasteten Scheibe 
anbringen muß, um den Lochrand lastfrei zu machen. 

Nach dem eingangs bewiesenen Integralsatz ver- 
schwindet aber das Integral der Störfunktion für einen 
Kreisweg nach Bild 2 solange der Radius des Integra- 
tionskreises als klein gegenüber den Abmessungen der 
Scheibe angesehen werden kann, weshalb 


$ (0% + 0), da = $ (0% + 0y)o da . 


Loch innerhalb des 
Integrationskreises 


Bas da= 1 DU 


. 2ald). 
lochfreie 
Scheibe 


Bild 3 


Null. Damit verschwindet aber das Integral GI. (1) 
auch für jede beliebige Randbelastung des Kreisloches. 
Dieser Beweis ließe sich auch bestätigen an Hand des 
streng berechenbaren Spannungszustandes einer un- 
endlich ausgedehnten Scheibe mit einem Kreisloch, 
dessen Rand beliebig belastet ist °). 


®) 8.2.B.A.u.L.Föppl: „Drang und Zwang“. Bd.I. 3. Aufl. 
8.286 ff., R. Oldenbourg, München 1941. 

®)S.z.B. K.Girkmann: „Flächentragwerke“, 8. 136 ff., 
Springer. Wien 1948. 


Bild 4 


Wenn das Loch vergrößert wird, bissein Rand mit dem 
Integrationsweg zusammenfällt (und damit zum Kreis- 
loch wird), dann gilt für den lastfreien Lochrand 

(% + oy)ı= 0: (6), 
worin o; die Tangentialspannung am Lochrand be- 
deutet. Damit folgt aus Gl. (5) 


Hada= Hgda (7), 


“% 8. 2. B. B.Baule: „Die Mathematik des Naturforschers und 
Ingenieurs“. Bd. VI S. 88, Hirzelverlag, Leipzig 1944. 
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wobei sich g(z, y) auf die ungelochte Scheibe bezieht. 
Wenn man beide Seiten dieser Gleichung durch 2 
dividiert, erhält man den Mittelwert der Randspan- 
nung des lastfreien Loches o+ in Identität mit der 
Spannungssumme Z im Mittelpunkt der Lochkontur 
der lochfreien Scheibe 


Dabei soll vorausgesetzt werden, daß die Abstände des 
Loches vom Rand der Scheibe großsind gegenüber den 
Abmessungen des Loches. 

Soweitstimmen diese Erkenntnisse mitdem Föppl- 
schen Mittelwertsatz überein, es wurden auch die 
gleichen Bezeichnungen dieser unter 1) zitierten Quelle 
gewählt. L.Föppl faßt seinen Mittelwertsatz aber 
viel weiter: 

„Der Mittelwert der Tangentialspannungen 0 
längs des Lochrandes einer unendlich ausgedehnten 
Scheibe stimmt mit dem Mittelwert der Normalspan- 
nungssumme g längs desselben Weges in der unge- 
störten lochfreien Scheibe überein. Einzige Voraus- 
setzung sollsein, daß die Störfunktion (x, y) keine 
Singularität außerhalb des Loches besitzt.‘‘ 

Die Gestalt des Loches ist danach keinerlei Einschrän- 
kung unterworfen. Diese Allgemeingültigkeit steht 
oder fällt mit der Aussage 


$ ot, y) = NT (9), 


wonach das Integral der Störfunktion über einen be- 
liebigen Lochrand verschwinden soll, wenn die Scheibe 
keine Singularitäten außerhalbdes Loches besitzt. Diese 
Störfunktion klingt mindestens mit dem Quadrat der 
Entfernung vom Loch ab, weshalb die Integration der 
Störfunktion über einen in großer Entfernung um das 
Loch geführten geschlossenen Integrationsweg den Wert 
Null ergibt. Da die Störfunktion zweifellos eine har- 
monische Funktion ist, sollauch für einen Integrations- 
weg, der näher um das Loch führt und schließlich auf 
seinen Rand zusammengezogen wird, das Integral 
der Störfunktion gemäß Gl. (9) verschwinden. 

Der Verfasser glaubt nun darin einen Irrtum insofern 
zu erkennen, als daß seines Erachtens für das all- 
gemeine geschlossene Bogenintegral einer reellen har- 
monischen Funktion keine physikalische Bedeutung 
in obigem Sinne besteht. 

Für eine komplexe analytische Funktion ist das 

B 


en! f(z) dz unter bestimmten Bedingungen 


vom Weg unabhängig. Nach dem Integralsatz von 
Cauchy verschwindet das geschlossene Linien- 
integral wenn keine Singularität eingeschlossen wird®): 


dfe)de=0 bzw. Hlplz, y) 
+iyla,y))(de +idyJ)=0. . (10a). 


Da Real- und Imaginärteil für sich verschwinden 
müssen, folgt 


gpde—yay=0; Hwde+pdy)=0 (1b), 


worin die Funktionen den Cauchy-Riemann- 
schen Beziehungen genügen müssen. In unserem Falle 
liegt nur die reelle harmonische Störfunktion (x, y) 
vor, undiihre Integration über einen geschlossenen Weg 
gemäß G]. (9) steht zu dem Linienintegral der Gl. (10) 
in keiner Beziehung. 

Es läßt sich leicht ein krasser Fall zeigen für den der 
Mittelwertsatz nicht gilt. 

In einer unendlichen Scheibe herrsche der allseitige 
Zug-Spannungszustand 0, = 0y = p. In diese Scheibe 
werde ein Loch geschnitten. Um die Lösung für den 
lastfreien Lochrand zu finden, müssen wir diesem 
gleichmäßigen Spannungszustand eine Störfunktion 
(0% + 0y) = P(x, y) überlagern, die hier von einem 


5) Z.B. R.Rothe: „Höhere Mathematik“, Bd. 2, S. 155 ff, 
u. 8.146ff., Teubner-Verlag Berlin 1942. 
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gleichmäßigen Druck —p auf den Lochrand herrührt. 
Das Loch habe die Gestalt von Bild 4, dann gilt für den 
äußeren in ausgezogener Linie gezeichneten Randteil 
%= — 04 bzw. og = — oy und eine Integration über 
diesen Randteil oder einen weiter außen angedeuteten 
Kreisweg ergibt Null. Für das Gebiet A einschließlich 
seinem gestrichelten Rand gilt aber = oy= —p. 

Die Störungen am Übergang dieser beiden Randgebiete 
können vernachlässigt werden. Das Wegintegral der 
Störfunktion p(z,y) über den gesamten Rand des 
Loches versch windet jetzt nicht mehr, sondern wird in 
vorliegendem Beispiel wegen der Länge des Weges um 
das Gebiet A einen beträchtlichen Wert annehmen. 

Das Verschwinden des Wegintegrals der Störfunktion 
%(x,y) über den Lochrand wäre aber Voraussetzung des 
Mittelwertsatzes. 

Wenn die tatsächliche Lochkontur größtenteilsin der 
Nähe eines Kreises verläuft, wird der Integrationswert 
der Störfunktion über diesen Lochrand nur wenig von 
der Integration über den Kreisweg abweichen, also an- 
nähernd verschwinden. Für diese Fälle stellt der 
Föpplsche Mittelwertsatz eine gute Näherung dar, 
wie spannungsoptische Versuche gezeigt haben®). 

Esläßtsich weiter zeigen, daß für jede Lochform eine 
unendliche Vielfalt zugeordneter Belastungen 
existiert, für die der Mittelwertsatz gilt. Diese Arbeiten 
des Verfassers sind noch nicht abgeschlossen. 


München. G. Sonntag. 


*e) L.Föppl: „Ein Beispiel zum Mittelwertsatz der ebenen 
Elastizitätstheorie‘‘. Z.angew. Math. Mech. 33 (1953), 8. 127—130. 


Numerische Lösung des Integrals 


Sflu)Itw(u) — w(g)? du 


l. Zusammenfassung. 
Integrale der Form 


ul Sa Hain) ze wa}? de ne 


erscheinen in der Theorie der laminaren Grenzschicht. 
Sie stellen Elemente von Lösungsmethoden dar, die die 
Stabilität der laminaren Grenzschicht bzw. den meist 
unerwünschten Umschlag in die turbulente Strömung 
zum Gegenstand haben!). Die zahlenmäßige Aus- 
wertung gestaltet sich oft schwierig. weil die Funk- 
tionen fu) und w(w) gewöhnlich nur für Gruppen 
diskreter reeller Werte bekannt sind, während die 
Lösung der komplexen Integrale Jeine gewisse Kennt- 
nis der analytischen Fortsetzung der Integranden vor- 
aussetzt. Gegenstand der vorliegenden Arbeit sind 
zwei Lösungsverfahren, die nur sehr wenig über die 
zahlenmäßigen Werte von f(w), w(u) und deren Ab- 
leitungen voraussetzen. Abgeleitet auf der Grundlage 
theoretischer Erwägungen, stellen sie doch eigentliche 
numerische Methoden dar, die überall dort geeignet 
sein dürften, wo die Funktionen f(w) und w(w) durch 
Verfahren der numerischen Analysis ermittelt wurden. 


2. Definition des Problems. 


Es soll angenommen werden, daß f(w) und w(u) die 
folgenden Eigenschaften besitzen: 


(a) Im Intervall a <u <b sind f(u) und w(u) ein- 
deutige, kontinuierliche und reelle Funktionen der 
reellen Veränderlichen «. 


(b) fw) >0 in (a,b). Sollte Su) das ee 
wechseln, wären Teilintegrationen in Teilintervallen 
von (a, b) vorzunehmen. 


4) Integrale vom Typus J treten in den folgenden theoretischen 
Abhandlungen auf: 


(a) Heisenberg W., „Über Stabilität und Turbulenz von Flüssig- 
keitsströmen‘‘, Annalen d. Physik, 24, 577, (1924). 
b) Lin C.C., „On the stability of two-dimensional parallel flows‘* 


Quarterly of Applied Mathematics, 3, Juli und Oktober 1945, 
(Jänner 1946). 
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(c) fw) und w(w) besitzen analytische Fortsetzungen 
in einen gewissen einfach zusammenhängenden Be- 
reich G, der das Geradensegment (a,b) in seinem Innern 
enthält. Diese Fortsetzungen sind jedoch nicht be- 
kannt. Bekannt sind nicht einmal analytische Aus- 
drücke für f(w) und w(w) im reellen Intervall (a,b). 
Die beiden Funktionen sind numerisch gegeben, für 
eine Gruppe äquidistanter reeller Werte von (a,b). 

(da<g<y<sb. 

(e) Der Integrationsweg C führe von einem Punkt A 
nach einem Punkt Y,in der Art und Weise, wie aus den 
Figuren (1) und (2) ersichtlich. 

Der Integrand von J besitzt einen Pol an der Stelle 
%= q. Wie unmittelbar ersichtlich, sind die beiden 
Konturen der Abbildungen (1) und (2) gleichwertig. 


u- Ebene 


u- Ebene 


Bild 1 


3.Lösungsmethode (I, Beseitigungder 
Singularität. 


Die hier entwickelte Methode ist auf dem Grundsatz 
aufgebaut, die Singularität des Integranden durch ana- 
lytische Umformungen zu beseitigen und das Integral J 
damit in eine Form überzuführen, die zahlenmäßige 
Auswertung nach konventionellen Verfahren (z.B. 
mit der Regel von Simpson) gestattet. 


Von der als bekannt vorausgesetzten Variation von 
f(u) und w(w) unmittelbar rechts und links von der 
Stelle = g können Näherungen für die ersten Koeffi- 
zienten der Ta ylorschen Reihenentwicklungen für 
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fu) und w(w) abgeleitet werden. Ist r eine sehr kleine 

positive Größe, und |u — q |< r, werden nur die aller- 

ersten Koeffizienten tatsächlich benötigt werden. 
uw)=Ahthu—- NMt+Ahu— gt. ), 
fW)=mtmwu—N)tmu— nt. @- 

Setzen wir (u — g) = t,erhalten wir 

Su) [ww — u = (UN wo + rnit mt?) 


worin die »; durch die A; und u; bestimmt sind. U. a. 
haben wir 2 


vo= Hol tw ()}° ale sale 
= — 100 (dw (+ Mltw (MP. (0). 


Die Reihenentwicklung (4) kann als Näherung der ana- 
lytischen Fortsetzung*des Integranden in eine kleine 
kreisformige Region um u= g aufgefaßt werden. 
Offenbar handelt es sich um die Laurentsche 
Entwicklung. Wir schreiben nun 


een (m) 


worin sich J, auf die Integration von Anach M bezieht 
(Abb. 2), J, auf die Integration von M nach N und J,. 
auf die Integration von N nach Y. J,ist leicht aus- 
gewertet. Man erhält 


1,0 
= —2nlr -inn +&2{1-(—1)% (1/n) 9,417” (8). 


Um nun auch J, und J, zu lösen, setzen wir 


Vf@)/ {w(u) — u(g)} 


= V{f}w (Du )+Iwg - - 9 
und somit 
 f® 
dr Hwi(g) 
Nay um, ; 
) nd er 
UM rg )- ra RO: 
g(u, = unse 
_V FL} wg) (u — D— tuta) — UT (11). 


w’(g) (u — g)-{w(u) — w(g)} 


Wie aus (11) und (12) ersichtlich, besitzt g(w, q) keine 
Singularität an der Stelle = g. Wir bilden das 
Quadrat von (9) und erhalten 


Stu) tw (u) — wid}? = <FO Hu (DI u — N? 
+ <U tfO}glu,g)/w’(Q> u — q) 
+ {g(u, 9}? (12). 


Um auch die Singularität des zweiten Summanden von 
(12) zu beseitigen, setzen wir 


gau,N)= gu) +hwg)-w—g.. . (13) 


mit dem Effekt, daß auch h(u, g) an der Stelle u=q 
nicht mehr singulär sein kann. Es ist 


h(g,q) = Og(u, q)[@u ander Stelleu= gq. . (14) 


Auf der Grundlage dieser Transformationen erhält 
man 


EN Rn A 


ae N 
(Mu —g) 


:9(9q) 1— 


gr „2/ ra) 
5 w’(q) 


-log|g—u| 


a 


ein) 


qa—r qa—r 
| Mn au + | {o(u, Q} du . (15) 


mn 
-—— 
— 


——— _...... 
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are ce AKIR Dea 
tw dPw—g|g+r w’(q) AL Pu, g) = <Ol tw} > (Du —d) 
r + pılq)log [u — gl 
2 Ua | + Palg) u 
a | mode + [rw ora ..a, Hug) 2 ge 
g+r g+r +. 
Summieren wir J), Ja und J,, fallen alle numerisch + <Ppn(Q ln — Dyila — , 0 +. (22). 


ungünstigen Ausdrücke heraus, und für den Grenz- 
übergang r — Oerhalten wir 


RE). 7, 


= (ia AI <Fd — SQ w’’(Dlw(Q> 
21/9: 2/7) N N 


DEREN 
RNIT LEN Ay — ga +1/g— a) 


Vu@} 


ne Er 


h(u, q) du 


+ [ (au, }* Zu — i fu kw) —w(g}? du (17). 


Mit der einzigen Ausnahme von h(q, q) gemäß (14) ent- 
hält kein Ausdruck von (17) Ableitungen von w(u) von 
höherem als dem zweiten Grad. Und diese einzige dritte 
Ableitung von w(u) kann leicht vermieden werden. 
Wir haben ja die Näherung 


g+e 
Sıa,)du=gats)—sa—-sd (IB). 


0 
vorausgesetzt, daßesehr kleinist. Alles was zur zahlen- 


mäßigen Lösung von (1) benötigt wird, ist demnach 
Kenntnis von: 


a) f(u) und w(h) für gewisse reelle Werte von win 
a, ’ 

b) Näherungen für die ersten Ableitungen von f(k) 
und w(u) für diese Stellen %;, 

e) Näherungen für w’’(w); w’’(u;) erscheint jedoch 
überall nur im Zähler und nicht im Nenner und ist 
demnach nicht mit besonders hoher Genauigkeit 
benötigt. 


4. Lösungsmethode (Il), Lösung durch 
analytische Fortsetzung. 

Unter gewissen Umständen könnte die folgende 
Lösungsmethode vorteilhaft sein. Wir nehmen eine 
Gruppe von Werten %, i=1,2,2,---, (n+]1) ge- 
geben an und setzen für einen dieser Werte = q. Ein 
Polynom 


F(u, q)= po(q) + Pıu — q) 
+91 — g)?+---pn(q) u — q)r (19) 


kann sodann aufgebaut werden nach den folgenden Be- 
stimmungen: 


wi) (wi — D’Kwlm) — wg): - 


i= 1,3,3,2.. —l), 
F (wi, = ( = 1), ... (n £ E 1), N (20), 
fdKw(dP? = Fi d 
= plN:.-t-ij 


F(u, qg) kann als Näherung für die analytische Fort- 
setzung von f(u) (u — g)?| {w(u) — w(g)}? betrachtet 
werden. Die durch Division mit (w — 9)? entstehende 
Funktion kann sodann als Näherung für die analytische 
Fortsetzung des Integranden von (1) betrachtet werden 
und durch gliedweise Integration erhält man 


J= P(y,)-P(ia)—pı)in . . (2l), 


Es ist ersichtlich, daß die einzige in der endgültigen 
Lösung von (l) benötigte wo von w(u) durch 
ihre erste gegeben ist. 


5.Lösung von (l) für große Werte vony. 

Wir wollen uns kurz der Frage zuwenden, wie (1) für 
große Werte von y zahlenmäßig gelöst werden kann. 
Dabei setzen wir der Einfachheit halber f(u) = 1 vor- 
aus. Die hier skizzierten Umformungen lassen sich 
leicht auf den allgemeineren Fallerweitern. Es soll also 


u 
fr du/tw(u) — w(g)} * bestimmt werden, wobei g, d, u 
u 


und u, der folgenden Abbildung entsprechen sollen. 
Zunächst wird angenommen, daß das Integral J für 
das Teilintervall (a, w) bereits nach einer derin den vor- 
angehenden Abschnitten behandelten Methoden aus- 
gewertet wurde. Wır setzen 
lim w(u) = 9 
ur”® 


Du Er Tore 29 


schreiben 


ie Ve Gi (24) 


und betrachten die asymptotische Entwicklung von 
w(u) in bezug auf d 


1 1 1 1: 
Wet ganz tgent . + (25). 
Setzen wir ferner 
w(u) = w(d + 1/2) = W(d, x) . (26) 
erhalten wir 
Wld,2)=n+t (ala + (gl2le?+--- . (27) 
4n = lim {on W(d, x)/0an . . . (28) 
>00 
und schließlich 
Us 
JS du kw(u) — w(g)}? 
u, 
1/(u, —d) 1 
1 
» (9 — w(q))* | zit 1784132? + .}dx (29), 
1/(u, —d) 


wobei die Koeffizienten r; von den Koeffizienten g; ge- 
funden werden, gemäß 

1/X1+(q1 a1! +92 @?]214 95 @®l3! ++ Hg weg)}22 

Sl en scno mo - (30). 

Es möge uns gestattet sein, die Brauchbarkeit der 

hier entwickelten Verfahren durch zwei Beispiele zu 


illustrieren. Dabei beschränken wir uns auf die Lö- 
sungsmethode (I). 


6. Beispiele. 

1. Beispiel. 

Nehmen wir an J, gegeben durch Gleichung (1), sei 
zahlenmäßig zu ermitteln, für fu) = 1 in (a,b), 
ww) =) win (a,b, y=1=b,a=1. 

J kann in diesem Fall exakt integriert werden und 
man erhält 

I = 2108 {1 -Yo/ d 21 —-Vg) —2ri al). 
Die für die numerische Auswertung entscheidenden 
Funktionen g und hergeben sich, andererseits, als 


wN)= 1a t/d; as d=URVd - - (82), 
32 


TE 
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"u, = - U2Vanda + 9°; 

h(q, q) = — (1B)IKAYd 
Die Lösung erscheint in der Gestalt. d 
J= 210g {1 I —AA a +21 +/D 

—2.log{1 +YdNF—2ri . (84). 

Es ist unmittelbar zu ersehen, daß (34) identisch ist mit 
(31). Größeres Interesse sollte jedoch der Beobachtung 
der äußerst günstigen numerischen Eigenschaften der 
Funktionen g(w,g) und h(u, qg) zugewendet werden. 
Beide Funktionen sind stetig, monoton und besitzen 
keinerlei Singularitäten. Das folgende, rein numerisch 
durchgerechnete Beispiel dürfte die grundlegenden 


Eigenschaften der hier entwickelten Lösungsverfahren 
jedoch noch deutlicher hervortreten lassen. 


2. Beispiel. 


Es sei g= 0.64, a= 0, b=y=1 und die Funk- 
tionen f(w) und w(w) seien zahlenmäßig gegeben durch 


u lee wu) 
0.0 0.0 0.0 
0.04 0.04 0.2 
0.08 0.08 0.282843 
0.12 0.12 0.346410 
0.16 0.16 0.4 
0.20 0.20 0.447214 
0.24 0.24 0.489898 
0.28 0.28 0.529150 
0.32 0.32 0.565685 
0.36 0.36 0.6 
0.40 0.40 0.632456 
0.44 0.44 0.663325 
0.48 0.48 0.692820 
0.52 0.52 0.721110 
0.56 0.56 0.748331 
0.60 0.60 0.774597 
q= 0.64 f(q) = 0.64 w(g) = 0.8 
0.68 0.68 0.824621 
0.72 0.72 0.848528 
0.76 0.76 0.871780 
0.80 0.80 0.894427 
0.84 0.84 0.916515 
0.88 0.88 0.938083 
0.92 0.92 0.959166 
0.96 0.96 0.979796 
1.0 1.0 1.0 


Außerdem seien die folgenden Ableitungen an der Stelle 
q gegeben: f’(g)= 1.0, w’(g) = 0.625000, w’’(g) = 
—0.488281. Es ergeben sich sodann gemäß der Be- 
stimmungsgleichungen für g und h die folgenden Zah- 
lenwerte (siehe folgende Tabelle). 


Durch Anwendung der Sim psonschen Regel erhält 
man weiter: 


1 
S {g(u, )}? du— 2.45410, 
0 


1 
S h(u, g) du= —0.256443 
0 


und schließlich, gemäß (17), 
J= S fu) wu) — w(g)}? du = — 7.522886 
c 
—3.84n.i\. 


Dem Leser wird vielleicht schon aufgefallen sein, daß 
die hier zahlenmäßig gegebenen Funktionen f(w) und 
w(w) identisch sind mit den analytisch exakten Funk- 
tionen 

Sw) 


vu und ww)=|u. 


Mitteilt N 


0.0 
0.04 
0.08 } 
0.12 ; 
0.16 s 
0.20 , 
0.24 , 
0.28 3 ’ 
0.32 2 
0.36 f 51 
0.40 ä 243 
0.44 i 233 
0.48 s 
0.52 1.52593 " | _2.32846 —0.216083 
0.56 1.51669 2.30035 —0.208625 
0.60 1.50806 2.27424 —0.201500 
gq= 0.64 1.5 2.25 —0.195500 
0.68 1.49242 2.22732 —0.189500 
0.72 1.48528 2.20606 —0.184000 
0.76 1.47853 2.18605 —0.175583 
0.80 1.47214 2.16720 —0.174125 
0.84 1.46606 2.14933 —0.169700 
0.88 1.46028 2.13242 —0.165500 
0.92 1.45476 2.11633 —0.161571 
0.96 1.44949 2.10102 —0.157843 
1.0 1.44444 2.08642 —0.154322 


Wir können nun vergleichsweise die exakte Integration 
vornehmen und erhalten: 


J= f uw — 0.8)? du 


= ln 
=1+4/9—2gl —Vg +00] 
q 
—bgrn.i 
— — 7.523354 — 3.834 -i, 


in zufriedenstellender Übereinstimmung mit dem Re- 
sultat der Näherungslösung. 


7.Schluß. 


Das Ziel der nun abgeschlossenen, sehr einfachen 
mathematischen Entwicklungen war, wie nochmals 
hervorgehoben werden möge, die zahlenmäßige Lösung 
von komplexen Integralen der Gestalt (1), auf der 
Grundlage numerisch gegebener Funktionen f(w) und 
w(u). Dies ist unzweifelhaft der in der praktischen 
Durchrechnung von Grenzschichtstabilitätsproblemen 
auftretende Fall. Mit Rücksicht darauf, daß die Na - 
vier-Stokesschen Gleichungen der inkompressib- 
len zähen Flüssigkeiten nur in Ausnahmefällen exakte 
Lösungen gestatten und es auch mit den Prandtl- 
schen Grenzschichtgleichungen nicht viel besser be- 
stelltist, kann die Funktion w(w) — die der Geschwin- 
digkeitsverteilung der Strömung entspricht — fast 
niemals anders als zahlenmäßig bestimmt werden. Das 
gleiche gilt sodann für f(w). In dieser Perspektive ge- 
sehen, darf vielleicht die Hoffnung ausgesprochen wer- 
den, daßsich die hier entwickelten Methoden praktisch 
brauchbar erweisen werden. 


Der Autor wünscht an dieser Stelle Herrn Professor 
Goertler vom Mathematischen Institut der Uni- 
versität Freiburg für die erhaltenen außerordentlich 
wertvollen Hinweise und Anregungen seinen besonderen 
Dank auszusprechen. 


Hawthorne (Los Angeles) 
California. Max Denglerh). 


1) Forschungsingenieur der Abteilung für Grenzschichtforschung 
von Northrop Aircraft Inc., Los Angeles, U.8.A. 
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BUCHBESPRECHUNGEN 


Jenaer Jahrbuch 1953. (Wissenschaft- 
liche Veröffentlichungen des Zeisswer- 
kes). Herausgegeben vom VEB Carl Zeiss Jena. 
269 S. m. zahlr. Abb. Jena 1953. Komissionsverlag 
VEB Gustav Fischer. Preis geb. 18.— DM. 

Auch das jetzt zum vierten Mal erscheinende Jenaer 
Jahrbuch bringt im wesentlichen Arbeiten derim VEB 
Carl Zeiss Jena beschäftigten Wissenschaftler. Die 
ersten drei Aufsätze sind der Geschichte des Zeiss- 
Werkes gewidmet; zwei würdigen die Verdienste von 
H. Hartig und A. König, früher verstorbener 
Mitarbeiter des Werkes, während der dritte die über 
hundertjährige Entwicklung des Zeiss-Mikroskop-Sta- 
tivbaues schildert. Die größere Zahl, nämlich sieben 
der Aufsätze, beschäftigen sich mit Fragen der reinen 
und der angewandten geometrischen Optik. Den Ab- 
schluß bildet ein theoretischer Aufsatz über die Dicken- 
schwingungen piözoelektrischer Platten und ein expe- 
rimenteller Beitrag über Dezimeterwellenantennen. 
Das Buch legt Zeugnis für die Pflege der Wissenschaft 
im Zeiss-Werk ab. 


Dresden. Willers. 


Prof. Dr. F,H. Müller, Das Relaxations- 
verhalten der Materie. (2. Marburger 
Diskussionstagung.) IV + 224 S.m. 112 Abb. Darm- 
stadt 1953. Verlag Dr. Dietrich Steinkopf. Preis 
brosch. 24.— DM. 

Nachdem auf der 1. Diskussionstagung 1950 der feste 
Zustand der hochpolymeren Substanzen behandelt 
worden war, befaßte sich die 2. Tagung mit den als 
Relaxationserscheinungen bezeichneten Zeiteffekten, 
die bei den in endlicher Zeit ablaufenden physikali- 
schen und chemischen Vorgängen auftreten, unter der 
Beschränkung, daß nur kleine Abweichungen vom 
Gleichgewicht auftreten. Es wird die thermodyna- 
mische Theorie der Relaxationserscheinungen auf der 
Grundlage der Thermodynamik der irreversiblen Pro- 
zesse dargestellt (Meixner) und hiermit die theore- 
tische Möglichkeit der Deutung mechanischer Relaxa- 
tionsspektren gezeigt (Staverman). Unter der 
Voraussetzung der Existenz mechanischer Relaxa- 
tionsspektren wird eine geschlossene, im wesentlichen 
aufeiner La place transformation beruhende mathe: 
matische Formulierung gegeben, die das Relaxations- 
spektrum aus dem Fließverhalten zu berechnen ge- 
stattet (Gr o ß). Ein statistisches Modell zur formalen 
Beschreibung von Relaxationsvorgängen wurde mit- 
geteilt (Ku bat). Verschiedene Vorträge befassen sich 
mit der experimentellen Feststellung (Schmieder 
u. W ol) diskreter Relaxationszeiten und ihrer theore- 
tischen Deutung sowie Auswertung unter Zugrunde- 
legung spezieller Fließformeln (Jenkel,Kneser). 
Ferner wurden u.a. para- und ferromagnetische Re- 
laxationsvorgänge und die dielektrische Dispersion 
(v.d. Marel; Würstlin; Brot, Magat u. 
Reinisch), neben Vorträgen über allgemeine Pro- 
bleme der spektralen Analyse und den Zusammenhang 
zwischen Platzwechselund Relaxation (H.F.Müller) 
behandelt. 


Aachen. F. Schultz-Grunow. 


Dr. Ing. Rudolf Richter (em. Prof. a. d. T.H. 
Karlsruhe) mit Beiträgen von Dr. Ing. Robert Brüder- 
link (o. Prof. a. d. T. H. Aachen). Elektrische 
Maschinen Bd. IL Synchronmaschi- 
nen u. Einankerumformer. 2. verb. Auf- 
lage XXVII 707 S. m. 525 Abb. Basel 1953. Verlag 
Birkhäuser. Preis geb. 46,80 Schw. Fr. 

Es zeugt besser als Worte vom grundsätzlichen Wert 
eines technischen Lehrbuches, wenn es, wie dieses, 
22 Jahre nach dem Erscheinen der ersten Auflage als 
zwar berichtigter und ergänzter, aberim wesentlichen 


gleichgebliebener photomechanischer Neudruck her- 
auskommen kann. Das diesem Buche seinerzeit von 
allen Seiten ausgesprochene Lob kann man auch heute 
noch auf die zweite vom Verlag Birkhäuser gut aus- 
gestattete Auflage übertragen. Das Kernstück ist 
nach einer kurzen Einführung über komplexe Zahlen, 
Ortskurven und Magnetfelder die Synchronmaschine 
inihrer Ausführung als Schenkelpol- und Vollpoltype. 
Neben dem Betriebsverhalten der Maschine als Einzel- 
maschine und in Zusammenarbeit miteinem Netz, den 
Ausgleichsvorgängen bei Kurzschluß und den Pende- 
lungserscheinungen, wird Anlassen, Synchronisieren 
und die prüffeldmäßige Untersuchung in ausführlicher 
und klarer Darstellung behandelt. Breiter Raum ist 
dem Entwurf eingeräumt, wobei dem heutigen Stande 
der Technik durch Ergänzen Rechnung getragen wird. 
Diese „Ergänzungen“ enthalten in gedrängter Form 
unter anderem Beiträge zur Kompoundierung von 
Synchronmaschinen, weisen auf die steigende Bedeu- 
tung des Synehronmotors hin, unterrichten über 
Haupt- und Hilfserregermaschinen, die Querfeldver- 
stärkermaschinen und über nichtstationäre Vorgänge 
unter Benützung derin Amerika gebräuchlichen Auf- 
teilung der Reaktanzen, sowie über dynamische Stabi- 
lität. An konstruktiven Neuerungen werden der 
Blechkettenläufer, die Kammpolkonstruktion und Ma- 
schinen mit Dauermagneten beschrieben. 

Der zweite Hauptabschnitt, der ebenfalls aus der 
1. Auflage übernommen wurde, gilt dem Einankerum- 
former. 

Es ist möglich, daß sich einige Leser die knappen 
Ergänzungen an einigen Stellen (z. B. über Wasser- 
stoffkühlung) ausführlicher wünschen, vielleicht sogar 
auf Kosten der Ausführungen über den heute durch 
die Gleichrichter stark verdrängten Einankerumformer. 
Man darf aber nicht vergessen, daß gerade diese tech- 
nisch interessante Maschinengattung einen guten Ein- 
blickin die elektrische Maschine erlaubt und deshalb 
für den Studierenden noch immer von Bedeutung ist. 
Außerdem gibt ein ausführliches Literaturverzeichnis 
über neue Arbeiten dem Leser leicht die Möglichkeit, 
sich ausführlicher zu unterrichten. 

Man kann überzeugt sein, daß die zweite Auflage 
auch weiterhin ihren guten Platz unter den Standard- 
werken der Elektrotechnik behaupten wird. 


Dresden. Pommer. 


Dr. Ing. Karl Kuhlmann (em. Prof. a. d. ETH 
Zürich. Theoretische Elektrotechnik 
Bd. II Grundzüge der Theorie’ elek- 
trischer Maschinen. (Lehr- und Hand- 
bücher der Ingenieurwissenschaften 13), 547 S. m. 
328 Abb. Basel 1951. Verlag Birkhäuser.. Preis geb. 
74,90 schw. Fr. 

Die als dritter Band eines Werkes über ‚‚Theore- 
tische Elektrotechnik“ erschienenen ‚‚Grundzüge der 
Theorie elektrischer Maschinen‘ sind in erweiterter 
Form aus den Vorlesungen des Verfassers an der Eid- 
genössischen Technischen Hochschule in Zürich her- 
vorgegangen. 

Das Buch bietetin einer vorzüglichen Ausstattung, 
für die dem Verlag Birkhäuser volle Anerkennung 
gebührt, mit seinen durch saubere Ausführung beson- 
ders ins Auge fallenden Zeichnungen in der Darstel- 
lung und Stoffauswahl eine eigene Arbeit, die von den 
gewohnten Lehrbüchern über elektrische Maschinen 
abweicht. Esliegt dem Verfasser offensichtlich weniger 
daran, Anleitung für den Entwurf und die Konstruk- 
tion zu geben, sondern vor allem in mathematischer 
Form das Verständnis für die grundsätzliche Wirkungs- 
weise elektrischer Maschinen zu vermitteln oder zu 
vertiefen. Es sind deshalb die dazu besonders geeigne- 
ten, zum Teil in großer Ausführlichkeit behandelten 
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Kapitel über die Erzeugung magnetischer Felder, die 
Induktivitätsberechnungen für Haupt- und Streu- 
induktivitäten und die induzierten elektromotorischen 
Kräfte wegen ihrer gründlichen Darstellung und der 
sehr guten Abbildungen lobend zu nennen. Ein Ab- 
schnitt mit der Überschrift Wicklungstechnik gibt 
Einblick in die verschiedenen Wicklungsarten, aller- 
dings kann man danach praktische Wieklungen wohl 
kaum entwerfen, ; 

Nach einem kurzen Abschnitt über magnetische 
Energie und das Drehmoment elektrischer Maschinen 
folgt ein Anhang von 152 Seiten, in dem besondere 
Vorgänge in den elektrischen Maschinen rechnerisch 
untersucht werden. Dort werden Kommutierungs- 
probleme, Fragen der Stromverdrängung und der 
Eisenverluste behandelt, außerdem wird das Rechnen 
mit symmetrischen Komponenten erläutert. Mehr 
betriebliche Fragen, wie Kennlinien, Kurzschlu Bver- 
halten und andere haben in dem Buch keine Auf- 
nahme gefunden. Die Art der Darstellung setzt, auch 
weil die vollständige Durchrechnung von einzelnen 
elektrischen Maschinen fehlt, ein gewisses Grundwissen 
über die elektrischen Maschinen bereits voraus und 
erfordert auch einige Bemühungen. 

Dem fortgeschrittenen Studierenden und dem wis- 
senden praktischen Ingenieur aber wird es eine Freude 
und eine gute Anregung für die eigene Arbeit sein, die 
Fülle der bisin die Einzelheiten theoretisch behandel- 
ten Probleme aus der Schau des Verfassers mitzusehen. 


Dresden. Pommer. 


Higher Transcendental Functions, 
Volume I. Based, in part, on notes left by Harry 
Bateman (Late Prof. at the California Institute of 
Technology) and compiled by the Staff ofthe Bateman 
Manuscript Projekt. XXVI +302S. New York/To- 
ronto/London 1953. Verlag Me Graw-Hill Book Com- 
pany, Inc. Preis geb. 56/6s. 

Der verstorbene Prof. Bateman plante ein um- 
fangreiches Werk über Theorie und Anwendung speziel- 
ler höherer Funktionen. Bei seiner umfangreichen 
Kenntnis des Schrifttums und seiner außerordentlichen 
Arbeitskraft würde er ein Werk geschaffen haben, das 
von der größten Bedeutung für die Anwendung dieser 
Funktionen in vielen Gebieten gewesen wäre. Als Vor- 
arbeiten dafür fand sich in seinem Nachlaß ein umfang- 
reicher Zettelkatalog. Sein umfassendes Programm 
erwies sich als nicht durchführbar und wurde von 
Erd&lyi, der unterstützt von W. Magnus, F. 
Oberhettinger, F. G. Triecomi die Heraus- 
gabe übernommen hatte, wesentlich gekürzt und auf 
die in den Anwendungen der Mathematik meist vor- 
kommenden Funktionen beschränkt. Dabei konnte in 
manchen Abschnitten nur teilweise von den Bate- 
manschen Aufzeichnungen Gebrauch gemacht wer- 
den, während sie für andere wesentlich benutzt wurden. 
Bei Funktionen von grundlegender Wichtigkeit und 
häufigem Vorkommen, wie es z. B. die Gamma-Funk- 
tionen sind, wurde mehr auf Einzelheiten eingegangen, 
wobei die Beweise allerdings auch da oft nur kurz ange- 
deutet wurden. In anderen Fällen ist die Darstellung 
mehr skizzenhaft oder es ist nur eine ausführliche Zu- 
sammenstellung der für den Gebrauch wichtigen For- 
meln gegeben. Dort wo ausführliche Darstellungen der 
früheren Ergebnisse vorliegen, wie etwa für die Bes- 
selschen Funktionen in dem Buch von Watson, 
konzentriert sich die Darstellung hauptsächlich auf die 
neueren Ergebnisse. Der vorliegende Band ist der erste 
des auf drei Bände berechneten Werkes über höhere 
transzendente Funktionen, dem noch zwei Bände mit 
Integraltafeln folgen sollen. Er behandelt die Gamma- 
Funktion, die hypergeometrische Funktion (Theorie 
und ausführliche Formelzusammenstellung), die Le - 
gendreschen Funktionen, die verallgemeinerten 
hypergeometrischen Reihen, ferner weitere Verallge- 
meinerungen der hypergeometrischen Funktionen, wie 


die E-Funktionen von Macrobert,die G-Funktin 
vonMeijer und die hypergeometrischen Funktionen 


mehrerer Variablen. Das letzte Kapitel ist der kon- 
fluenten hypergeometrischen Funktion gewidmet. Das 
‚Buch bringt trotz der erwähnten Kürzungen auch jetzt 
noch einen außerordentlich umfangreichen Stoff. Man 
muß den Bearbeitern für die auf die Herausgabe ver- 
wandte Mühe dankbar sein und darf nur hoffen, daß die 
fehlenden Bände nicht allzulange aufsich warten lassen. 


Dresden. Willers. 
> I 


Proceedings of the Eastern Joint Com- 
puter Conference. 125 S. m. zahlr. Abb. New 
York 1954. The Institute of Radio Engineers, Inc. 
Preis geh. 3 $. 

Auf der drei Tage dauernden Konferenz wurde in 
rund 25 Vorträgen, deren jeder über umfangreiche 
Untersuchungen berichtet, eine Übersicht über die 
Anwendung der heute in den USA arbeitenden, etwa 
50 größeren Ziffernrechenautomaten und die mit 
ihnen gewonnenen Erfahrungen, insbesondere was ihre 
Zuverlässigkeit betrifft, gegeben. So behandeln die 
ersten Vorträge die Verwendung auf dem Gebiet der 
Lebensversicherung, des Luftverkehrs, der Wetter- 
vorhersage usw. Andere geben die Erfahrung wieder, 
die man z. B. mit der REAC, der SEAG usw. oder 
auch mit Einzelteilen, wie mit den Magnetkernspei- 
chern, den elektrostatischen Speichern usw. gemacht 
hat. Der allgemeine Eindruck, den man aus diesen 
Vorträgen gewinnt, ist der, daßin den letzten Jahren 
außerordentliche Fortschritte gemacht wurden. Aus 
der Hauptentwicklungslinie und aus den einzelnen 
noch vorhandenen Lücken wird auf die wünschens- 
werte Weiterentwicklung geschlossen. Insbesondere 
sollte weniger Gewicht auf den Bau großer, allgemein 
brauchbarer Maschinen, als auf den von Maschinen 
für spezielle Zwecke gelegt werden, z. B. sollte man 
Wert auf den Bau von Maschinen mit Höchstge- 
schwindigkeit legen, wie sie für die Verwendung zur 
Wetterprognose oder für den Einbau in Regelungsvor- 
richtungen in Frage kommen. 


Dresden. Willers. 

Dr. Ing. R. Richter (Prof. em. a. d. T. H. Karls- 
ruhe), Elektrische Maschinen, Band III: Die 
Transformatoren. Zweite, verb. Aufl. 321 S. m. 
230 Abb. Basel/Stuttgart (1954). -Verlag Birkhäuser. 
Preis geb. 26 Fr. (26.—DM). 

Auch die 2. Auflage des Bandes 3 ‚Transformatoren “ 
des längst anerkannten dreibändigen Standardwerkes 
von Richter, wird von zahlreichen Ingenieuren der 
Praxis und vor allen Dingen auch vom Ingenieurnach- 
wuchs, den Studententen, hoch geschätzt werden. Die 
Neuauflage ist ein fotomechanischer Neudruck; durch 
die dem Buch vorangestellten Ergänzungskapitel und 
das umfangreiche zusätzliche Literaturverzeichnis wird 
trotzdem der neueste Stand der technischen Entwick- 
lung berücksichtigt. Wertvoll an dem Richterschen 
Werk ist die bekannte Gründlichkeit des Verfassers, 
mit der alle für die Tranformatoren wichtigen Teilfragen 
behandelt werden. 

Für das Kapitel K: ‚Beanspruchung und Wahl der 
Isolierung‘ möchte der Berichter, für den dieses Kapitel 
von besonderem Interesse ist, einen Hinweis machen, 
der bei der Beurteilung der Berechnungen von Trans- 
formatoren beachtet werden muß. £ 

Im Kapitel 3: „Zusammengesetzte Isolierstoffe‘‘ wird 
die Bemessung der Isolierung zwischen den Wicklungen 
eines Transformators behandelt. Die Zahlentafel 12 
gibt die auftretenden Beanspruchungen E;.der dünnen 
und Tafe] 13 die der dicken Isolierzylinder an, wenn 
das Öl leitend wird und der feste Isolierstoff die gesamte 
Beanspruchung auszuhalten hat. Die Zahlenwerte sind 
nach Gleichung 222a und b unter der Annahme eines 
nahezu homogenen Feldes berechnet worden. In Wirk- 
lichkeit werden jedoch diese Isolierzylinder dann, wenn 


Bun 


unse 
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das Ölleitend wird, d.h. teilweise durchschlägt, durch 
das Feld Spitze- Spitze d.h.durch ein stark inhomogenes 


Feld beansprucht, für das die in den Werkstofftafeln 


enthaltenen Festigkeitszahlen nicht mehr gelten. Es be- 
steht also dann nichtnur für die dünnen, sondern auch 
für die dieken Zylinder bei Teildurchschlag des Öles die 
Gefahr des Gesamtdurchschlages oder zumindest ernster 
Beschädigung auch bei kurzer Prüfdauer. 


Dresden. Obenaus. 


L. Howarth, H.B. Squire and C.N.H.Lock, Mo- 
dern Developments in Fluid Dyna- 
mies. High Speed Flow. (The Oxford 
Engineering. Science Series.) 2 Bände. XVI-+ VIII 
+ 8758. m. 288 Abb. u. 26 Tafeln. Oxford 1953. 
Clarendon Press. Preis geb. 84.—s. 

Die mit der Entwicklung des Hochgeschwindigkeits- 
flugzeugbaus, der Raketenflugtechnik sowie der Bal- 
listik in den letzten Jahrzehnten aufgetretenen viel- 
fältigen Probleme erforderten eine intensive For- 
schungsarbeit auf dem Gebiet der Dynamik kom- 
pressibler Strömungen. Die bisherigen allerseits unter- 
nommenen Anstrengungen in dieser Richtung haben 
eine Fülle von theoretischen und experimentellen Er- 
fahrungstatsachen erbracht. Das vorliegende Werk 
hat sich eine systematische Zusammenfassung dieses 
Materials zum Ziel gesetzt. Es bildet damit entspre- 
chend dem Wunsch der Herausgeber eine Ergänzung 
zu demim Jahr 1938 unter der Redaktion von 8. Gold- 
stein abgefaßten zweibändigen Handbuch ‚Modern 
Developments in Fluid Dynamics‘, in dem vornehm- 
lich die inkompressiblen Flüssigkeiten behandelt 
werden. 

Das Gelingen, das dem Werk in jeder Hinsicht be- 
scheinigt werden kann, ist wohl nicht zuletzt darin be- 
gründet, daß die Abfassung der 11 Kapitel, in die der 
umfangreiche Stoff gegliedert ist, jeweils einzelnen 
Autoren — Mitarbeitern des britischen Aeronautical 
Research Council — übertragen wurde, und zwar: 
Introduetion — L. Howarth ; The Equations of Flow 
in Gases —L. Howarth ; The Method of Characteristics 
— R.E. Meyer; Shock Waves —C.R. Illingworth; 
Blast Waves — G. J. Kynch; Some Exact Solutions 
of the Equations of Steady Homentropic Flow ofan 
Inviscid Gas — W. G. Bickley; One-dimensional 
Flow — 0.A. Saunders; The Hodograph Trans- 
formation — M. J. Lighthill; Approximate Methods 
— G. N. Ward; Unsteady Motion — G. Temple ; Boun- 
dary Layers — A.D. Young; - Experimental Methods 
— D. W. Holder, D.C. Macphall, J. S. Thompson; 
Flow Past Aerofoils and CGylinders — W.A. Mair, 
J. A. Beavan; Flow Past Bodies of Revolution — 
W.F. Cope; Heat Transfer — H.B. Squire. Da- 
mit konnte erreicht werden, daß jedes Spezialgebiet 
fachmännisch bearbeitet ist und alle wichtigen neuen 
Ergebnisse die ihnen gebührende Berücksichtigung 
finden, sei es in dem behandelten Stoff, sei esin den 
überaus zahlreichen Literaturhinweisen. Dabei ist 
noch die der Redaktion zu dankende Wahrung der 
einheitlichen Darstellungsweise hervorzuheben. Daß 
in den Literaturangaben vielfach wichtige, in der 
UdSSR durchgeführte Forschungsarbeiten nicht er- 
wähnt werden, dürfte wohl in erster Linie auf der 
schweren Zugänglichkeit dieser Arbeiten beruhen. 

Wenngleich auch häufig auf die in dem oben ge- 
nannten ersten Werk enthaltenen Ergebnisse Bezug 
genommen wird, so werden doch, um auch hier eine 
in sich geschlossene Darstellung zu erreichen, alle die 
Dynamik kompressibler Flüssigkeiten betreffenden 
allgemeinen theoretischen Grundlagen vollständig be- 
handelt. Das Werk ist daher sowohl als Lehrbuch wie 
auch als ein dem neuesten Stand entsprechendes 
Handbuch für den Fachman auf dem Gebiet der Aero- 
dynamik wärmstens zu empfehlen, letzteres besonders 
im Hinblick auf die ausführliche Behandlung der 
experimentellen Methoden und Ergebnisse. 


Dresden. RK. Krienes. 
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N. I. Muskhelishvili (Math. Institut Tiflis), Singu- 
lar integral equations. 2. Aufl. (Aus dem Russi- 
schen übersetzt von J.R.M. Radok.) 447 S. m. 16 Abb. 
Grooningen (Holl.) 1953. Verlag P. Nordhoff. N.V. 
Preis geb. 28,50 fl. 

Bei der Behandlung verschiedener Randwertauf- 
gaben der mathematischen Physik, vor allem der 
Potentialtheorie und ihrer Anwendungen auf Elastizi- 
tätstheorie und Hydrodynamik, stößt man auf 
singuläre lineare Integralgleichungen, deren Integral- 
glied als Cauchyscher Hauptwert zu verstehen ist und 
die sich also der Fredholmschen Theorie nicht mehr 
einfügen. Es ist das große Verdienst des Verf., alle 
diesbezüglichen Einzeluntersuchungen gesammelt 
und einer Theorie der (komplexen) Integralgleichung 
vom ‚„Cauchyschen Typ“ 


il K d 
A re m 
L 


(4A, K, f bekannt, p gesucht, Z der Integrationsweg in 
der komplexen Ebene) eingeordnet zu haben, die 
man dem Verf. und seinen Mitarbeitern (insbesondere 
I. N. Vekua) verdankt und die, wenigstens in dem 
hier vorliegenden Falle von Linienintegralen, eine ge- 
schlossene Darstellung gestattet. Dieser Theorie ist 
das vorliegende Werk gewidmet. 

Verf. verzichtet bewußt auf allzu große Allgemein- 
heit der Voraussetzung und begnügtsich mit gewissen 
Einschränkungen für gesuchte und gegebene Funk- 
tionen (Hölderbedingung), was zu abschließenden Er- 
gebnissen führt und dabei den Anforderungen der 
praktischen Anwendungen durchaus genügt. Im 
ersten Teil findet sich ein ausführliches und sehr sorg- 
fältiges Studium von Integralen der in (*) auftreten- 
den Art. Esfolgt dann die Auflösungstheorie, die auf 
der Zurückführung auf eine Fredholmsche Gleichung 
unter Verwendung geeigneter Integraloperatoren be- 
ruht und zu Resultaten führt, die den. Fredholmschen 
Sätzen verwandt sind. Dabei kann Z aus einer end- 
lichen Anzahl geschlossener Kurven oder auch, was 
die Betrachtungen wesentlich kompliziert, aus einer 
endlichen Anzahl von Kurvenstücken bestehen. An- 
wendungen bilden eine Reihe entsprechender Rand- 
wertaufgaben der Funktionentheorie, Potentialtheo- 
rie, Elastizitätstheorie und Hydromechanik für von 
geschlossenen Kurven berandete Gebiete oder für die 
längs Kurvenstücken aufgeschlitzte' Ebene. So fin- 
den u. a. auch das in der Theorie der dünnen Profile 
auftretende Umkebrproblem für Cauchysche Integrale 
und, nach einiger Umformung, die Prandtlsche Trag- 
flügelgleiehung innerhalb dieser Betrachtungen ihren 
Platz. Den Abschluß des Buches bilden entsprechende 
Untersuchungen für Systeme von zu (*) analogen 
Integralgleichungen 

Das Buch ist für jeden lesbar, der die Grundlagen 


‚der Funktionentheorie und die Lehre der Fredholm- 


schen Integralgleichungen kennt. Eine Übersetzung 
dieses ausgezeichneten Werkesin die deutsche Sprache 
wird vom Akademie-Verlag, Berlin, vorbereitet. 


Dresden. K. Maruhn. 


Dr. H. Blasius (Oberstudienrat i. R. der Ingenieur- 
schule Hamburg), Höhere Mathematik. Mathe- 
matische Grundlagen vom technischen Standpunkt. 
294 S. m. 168 Fig. und 78 Aufgaben. Hamburg 1954. 
Boysen & Maasch Verlag. Preis kart. 8, — DM, geb. 
10,— DM. 

Der Grundgedanke im methodischen Aufbau des 
vorliegenden Werkes ist die Feststellung, daß ab- 
strakte Begriffe und Sätze erst dann vom Leser richtig 
erfaßt werden können, wenn diese sich aus praktischen 
Aufgaben heraus sozusagen aufdrängen. Der Verf. 
will „zeigen, wie die Methoden zwecks Lösung von 
Aufgaben entwickelt sind‘ und daß, ‚die Begriffe aus 
dem Beispiel erwachsen müssen‘. Nach diesem Prin- 
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zip werden, ausgehend von Beispielen aus Geometrie 
und Mechanik, die im Untertitel genannten Gebiete 
in Angriff genommen. Der vorgesehene Stoff gliedert 
sich in 10 Gruppen, innerhalb deren von geeigneten 
Beispielen zum gewünschten mathematischen Gegen- 
stand aufgestiegen wird (z. B. von Geschwindigkeit 
und Steigung zur Differentialrechnung, von Fläche 
und Schwerpunkt zur Integralrechnung usw.). 

Daß die schließlich erarbeiteten Definitionen im 
allgemeinen der Anschauung entstammen und daß 
bei vielen wichtigen Sätzen die Beweise durch Plau- 
sibilitätsbetrachtungen ersetzt sind, ist an sich im 
Hinblik auf den Leserkreis (wohl in erster Linie Stu- 
dierende an Fachschulen) durchaus kein Unglück. 
Es müßte aber auf — nun einmal notwendige — Lük- 
ken ausdrücklich hingewiesen werden, um den kritisch 
eingestellten Leser einen kleinen Fingerzeig zu geben, 
was nun eigentlich ‚höhere Mathematik“ ist. 

Die geschickt ausgewählten Beispiele und Aufgaben 
werden gar manchem als Nachschlagewerk oder auch 
bei der Vorbereitung auf Prüfungen gute Dienste 
leisten. 


Dresden. K. Maruhn. 


W.Blaschke (Ord. Prof. a. d. Univ. Hamburg), 
Analytische Geometrie (Lehrbücher und Mono- 
graphien aus dem Gebiete der exakten Wissenschaften. 
Mathematische Reihe, Band 16). 2. Aufl. 190 S. m. 
67 Abb. Basel/Stuttgart 1954. Verlag Birkhäuser. 
Preis geb. 19,60 DM (19,60 Fr.), brosch. 16,65 DM 
(16,65 Fr.). 

Zu der 1. Auflage (Wolffenbütteler Notdrucke 1946 
Besprechung ds. Zeitschr. Bd.29 (1949), S. 62/63) fügte 
Verf.ein Kapitel über Waben und Gruppen hinzu. Er 
führt uns dort an Hand der von Reidemeister und 
Thomson entwickelten Wabentheorie in den Geist der 
modernen Axiomatikein. 


Halle. Ott-Heinrich Keller. 


W. Blaschke (Ord. Prof. a.d. Univ. Hamburg), Pro- 
jektive Geometrie (Lehrbücher und Monogra- 
phien aus dem Gebiete der exakten Wissenschaften. 
Mathematische Reihe Band 17). 3. erw. Aufl. 197 S. 
m. 71 Abb. Basel/Stuttgart 1954. Verlag Birkhäuser. 
Preis geb. 19,60 DM (19,60 Fr.), bresch. 16,65 DM, 
(16,65) Fr.). 

3. Auflage des überall geschätzten Bändchens aus 
der Reihe der Wolffenbütteler Notdrucke, (Bespre- 
cehungs. ds. Zeitschr. Bd. 28 (1948), S. 250/51). 


Halle. Ott-Heinrich Keller, 


A. Erdelyi, W. Magnus, F. Oberhettinger, F. G. Tri- 
comi (Staff of the Bateman Manuscript Projekt). 
Tables of Integral Transforms. Vol. I. (Based, 


in part, on notes left by Harry Batemann, Late Prof.- 


at the California Institute of Technology). XX + 
391 S. New York/Toronto/London 1954. Me Graw-Hill 
Book Comp., Inc. Preis 56/6 s. 

Dem ausdem Nachlaß von H. Bate man herausgege- 
benen Werk über höhere transzendente Funktionen, 
dessen erster Band auf Seite 476 angezeigt wurde, 
sollteein zweibändiges Tafelwerk von dem gleichen 
Herausgeberkreis unter Benutzung des von Bateman 
bereitgestellten Materials angeschlossen werden, von 
dem hier der erste Band vorliegt. Dieses Tafelwerk ist 
für den Gebrauch von Mathematikern, Physikern und 
Ingenieuren bestimmt. Es sollte die in der reinen und 
angewandten Mathematik häufig gebrauchten Inte- 
graltransformationen und sonstige höhere Integrale 
enthalten. Im ersten Bande finden sich in Form von 
Tabellendie Sinus-,Cosinus- und Exponential-Fourier- 
Transformierten, die Laplace-Transformijerten und die 
inversen Laplace-Transformierten, die Mellin-Transfor- 
mierten und die inversen Mellin-Transformierten. Die 


Anordnung der Tafeln ist 
nen Formeln zuerst die elementaren 
men (rationale, algebraische Funktionen, 
die Potenzen mit beliebigen nicht notwendi 
len Exponenten enthalten, Exponentialfunktion: 
Logarithmen, trigonometrische und hyperbolis 
Funktionen mitihren Umkehrfunktionen). Weiter fol- 
gen die höheren Funktionen (Orthogonal-Polynome, 
Gamma-Funktionen, Fehlerfunktionen, Besselsche 
und damit verwandte Funktionen, parabolische Zylin- 
derfunktionen, hypergeometrische Funktionen und 
ihre Verallgemeinerungen, elliptische Funktionen 
und einzelne sonstige Funktionen). In verschiedene 
Gruppen gehörige Funktionen sind meist mehrfach ww 
aufgeführt. Zusammen mit dem hoffentlich bald er- % 
scheinenden zweiten Band ist so ein außerordentlich a 
umfangreiches, übersichtlich geordnetes Tafelwerk a 
entstanden, das weite Verbreitung finden wird und 

für dessen sorgfältige, mühevolle Bearbeitung alle 

Benutzer den Herausgebern dankbar sein werden. 


Dresden. Willers. 


D. N. de G. Allen, Relaxation Methods. 
IX + 257 S. mit 133 Abb. New York/Toronto/London 
1954. McGraw-Hill Book Comp., Inc. Preisgeb. 53/6 s. 

Der Verfasser, ein langjähriger Mitarbeiter von 
Southwell, gibt hier eine ausführliche Erläuterung 
der Technik des Relaxationsverfahrens unter Be- 
nutzung vieler Beispiele. Dabei behandelt er nur 
einfache Probleme, gibt ferner im allgemeinen nicht 
die Ableitungen der Gleichungen, deren Lösung er 
durchführt und die meistens der Elastizitätstheorie, 
der Hydrodynamik und der Theorie der Elektrizität 
und des Magnetismus entnommensind,geht auch kaum 
auf die Theorie des Verfahrens ein, z. B. nicht auf die 
Fragen, die mit dem Ersatz der Differentialgleichun- 
gen durch Differenzengleichungen zusammenhängen, 
auf die Wirkung der Abrundungsfehler usw. Der Stoff 
istin 16 Kapitel aufgegliedert: Lineare algebraische 
Gleichungen, Fachwerkprobleme, gewöhnliche Dif- 
ferentialgleichungen, Laplace’sche und Poisson’sche 
Gleichungen, Grenzbedingungen, Quasi-ebene Poten- 
tialgleichungen, Biharmonische Gleichungen, simul- 
tane Differentialgleichungen, Dreiecknetze, Eigen- 
wertprobleme bei Matrizen und Differentialgleichun- 
gen, mehrfach zusammenhängende Bereiche, Bereiche 
mit noch festzulegenden Begrenzungen, nichtellip- 
tische Differentialgleichungen und schließlich drei- 
dimensionale Probleme. Besonders betont der Ver- 
fasser:einmal, daß man mit dem Verfahren nur dann 
erfolgreich arbeiten kann, wenn man es sich durch 
systematische Übung angeeignet hat, wenn man also 
etwa die am Ende der einzelnen Kapitel gegebenen 
Aufgaben der Reihe nach wirklich durchgerechnet hat, 
und zum anderen, daß das Relaxationsverfahren kein 
starres Rechenschema ist, sondern daß durch sinn- 
volle Verwendung der Punkt- oder der Blockrelaxa- 
tion, der Über- oder der Unterrelaxation sowie durch 
Anwendung des sogenannten multiplizierenden Fak- 
tors eine Mannigfaltigkeit in das Verfahren kommt, 
die zur möglichsten Beschleunigung der Konvergenz 
benutzt werden kann. 


Dresden. Willers. 


CGolloque sur les fonctions de plusieurs var- 
riables.Tenu & Bruxelles du 1lau 14 mars 1953.1598. 
Liöge, Paris 1953. Verlag Georges Thone, Masson & 
Cie. Preis geh. 1.800 francs francais, 250 francs belges. 


Table des matieres 


Severi, F., Quelques probl&mes se rapportant aux 
fonctions analytiques de plusieurs variables 

Lelong, P., Fonctions plurisousharmonigues; me- 
sure de Radon associ6es. Application aux fonctions 
analytiques. 
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Cartan, H., Vari6t6s analytiques complexes et co- 
homologie. 

Serre, J.-P., Quelques probl&mes globaux relatifs 
aux vari6t6s de Stein. 

Roquette, P., L’Arithmötique des fonctions 


' aböliennes. 


Behnke, H., Gön6ralisation du thöor&me de Runge 
pour des fonctions multiformes de variables com- 
plexes. 

Stein, K., Analytische Projektion komplexer Man- 
nigfaltigkeiten. 

Martinelli ‚E .,‚Sur l’extension des th&or&mes de 
Cauchy aux fonctions de plusieurs variables com- 
plexes. 

Saxer, W., Sur les domaines de normalit6 des fonc- 
tions möromorphes de plusieurs variables. 

Bergman, S., Kernel function and extended Clas- 
sesin the theory of functions of complex variables. 


Dresden. K. Maruhn. 


Premier Colloque sur les öquations aux 
d6öriv6es partielles. (Tenus & Louvain du 17 au 
19 d&cembre 1953). Centre Belge de recherches mathe- 
matiques. 123 S. Liege, Paris 1954. Masson & Cie, 
G. Thone. Preis geh. 200 francs belges. 1400 francs 
frangais. 

Table des matieres 

Lichnerowicz, A., Equations de Laplace et 
espaces harmoniques. 

Four&s, Y., Rösolution de problöme de Cauchy 
pour les öquations hyperboliques du second ordre non 
lineaires. 

Delsarte, J., Sur certains syst&mes d’dquations 
aux derivöes partielles & une seule fonctioninconnue,et 
sur une gön6ralisation de la thöorie des fonctions de 
Bessel et des fonetions hyperg&omötriques. 

Doetsch, G., L’application de la transformation 
bidimensionnelle de Laplace dans la thöorie des öqua- 
tions aux d6riv6es partielles. 

Lepage, Th., Equations du second ordre et trans- 
formations symplectiques. 

Gillis, P., Sur certaines classes d’&quations aux d6ri- 
vöes partielles du second ordre, non linöaires. 

Sauer, R., Remarques gdomötriques sur les 6qua- 
tions aux dörivses partielles du second ordre quasili- 
nöaires et homogönes. 


Dresden. K. Maruhn. 


N.I.Muskhelishvili, Some basic problems 
ofthe mathematical theory ofelasti- 
eity. 3. verb. und vermehrte Auflage. (Aus dem 
Russischen übersetzt von J. R.M. Radok.) XXXI + 
704 S. m. 66 Abb. Groonigen (Holl.) 1953. Verlag 
P. Noordhoff Ltd. Preis geb. 38,— fl. 

Das Erscheinen einer englischen Übersetzung der 
dritten Auflage der preisgekrönten Monographie von 
Muskhelishvili über mathematische Elastizitäts- 
theorie wird sicher in den englisch sprechenden Län- 
dern mit Begeisterung begrüßt werden. Vor etwa 
20 Jahren wurde eine erste Auflage des Werkes in 
russischer Sprache veröffentlicht, welche jedoch aus 
sprachlichen Schwierigkeiten weiten Kreisen unbe- 
kannt geblieben war, so daß die allgemeinen Lösungs- 
methoden für zweidimensionale Probleme der Elasto- 
statik, wie sie durch Muskhelishvili entwickelt worden 
sind, bisher noch nicht den gebührenden Platz in der 
Weltliteratur eingenommen hatten. Inhaltlich um- 
faßt das Buch folgende Teilgebiete: Grundlegende Be- 
ziehungen der Mechanik elastischer Körper (Span- 
nungsanalyse, Formänderungsanalyse, Grundgesetze 
der Elastizitätstheorie).. Allgemeine Formeln der 
ebenen Elastizitätstheorie (Grundgleichungen der 
ebenen Elastizitätstheorie, die Spannungsfunktion 
und die komplexe Darstellung der allgemeinen Lösung 
der ebenen Elastizitätsgleichungen, mehrdeutige Ver- 


‚s—. 
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schiebungen und Wärmespannungen, Transformation 
der Grundbeziehungen für konforme Abbildung). 
Lösung einiger Probleme der ebenen Elastizitäts- 
theorie mit Hilfe von Potenzreihen (Fourier-Reihe, 
Lösungen für Kreisgebiete, Kreisring, Anwendung der 
konformen Abbildung). Cauchy-Integrale (Grund- 
eigenschaften von Cauchy-Integralen, Randwerte 
holomorpher Funktionen). Anwendung der Cauchy- 
Integrale zur Lösung von Randwertproblemen der 
ebenen Elastizitätstheorie (allgemeine Lösung der 
Grundprobleme für einfach zusammenhängende Be- 
reiche, Lösung der Grundprobleme für Gebiete, welche 
durch rationale Funktionen auf einen Kreis abgebildet 
werden können, Erweiterung auf eine Näherungs- 
lösungfür Gebiete allgemeiner Art, Lösung der Grund- 
probleme für die Halbebene und für halb-unendliche 
Gebiete, allgemeine Lösungsmethoden für Randwert- 
probleme). Lösung des Problems der ebenen Elastizi- 
tätstheorie mittels Zurückführung auf das Problem 
der linearen Verwandtschaft der Randwerte (Lösung 
von Grundproblemen für die Halbebene und für die 
Ebene mit geradlinigen Ergänzungen, Lösungen von 
Randproblemen für Gebiete, welche durch Kreise be- 


grenzt sind, und für die unendliche Ebene, welche - 


durch Kreisbögen geschnitten wird, Lösung der Rand- 
probleme für Gebiete, die sich durch rationale Funk- 
tionen auf dem Kreis abbilden lassen). Zug, Torsion 
und Biegung homogener und zusammengesetzter 
Stäbe (Torsion und Biegung homogener Stäbe, Tor- 
sion von Stäben aus verschiedenem Material, Zug und 
Biegung von Stäben aus verschiedenem Material mit 
gemeinsamer Poisonscher Konstante, Zug und Bie- 
gung bei verschiedenen Poisonschen Konstanten). 
Anhang; Begriff des Tensors, Bestimmung von Funk- 
tionen aus ihren Differentialen in mehrfach zu- 
sammenhängenden Gebieten, Bestimmung . einer 
Funktion einer komplexen Variablen aus dem Real- 
teil, unbestimmte Integrale holomorpher Funktionen. 
Literaturverzeichnis, Sachverzeichnis. Das Werk 
bringt eine weittragende Methode zur Integration der 
Grundgleichungen der ebenen Elastizitätstheorie mit 
Hilfe komplexer Funktionen. Die auf anderem Wege 
nur mit großen Schwierigkeiten durchführbare Lö- 
sung ebener Spannungsaufgaben wird durch den Ver- 
fasser weitgehend vereinfacht und für die praktische 
Anwendung zugängig gemacht. Das Buch kann nicht 
nur allen Mathematikern, die sich mit der mathe- 
matischen Elastizitätstheorie beschäftigen möchten, 
sondern auch jedem mathematisch vorgebildeten 
Ingenieur, der sich für höhere Festigkeitslehre in- 
teressiert, wärmstens empfohlen werden. 


Dresden. H. Neuber. 


Dr.-Ing. M. Schuler (Prof. em. der Universität Göt- 
tingen, Einführung in die Mechanik 
Teill. Mechanik des Massenpunktes. 
2028. m. 114Abb. Teil II. Mechanik der 
Punktsysteme. 2108. m. 166 Abb. Wolfen- 
büttel 1950/51. Verlag Benno Kracke. Preis geb. 
je 12,— DM. 

Verfasser veröffentlicht mit dem vorliegenden Werk 
eine Einführung in die Mechanik, die aus einer Reihe 
von Vorlesungen entstanden ist, welche von ihm in 
den Jahren 1926 —1946 an der Universität Göttingen 
gehalten wurden. Es wird die Kenntnis der Grund- 
lagen der Differential- und Integralrechnung, sowie 
der elementaren Vektorrechnung vorausgesetzt. Der 
Inhalt ist allgemein verständlich geschrieben. 

Der erste Teil gliedert sich in folgende Abschnitte: 
Der Newtonsche Satz und seine Begriffe; die gerad- 
linige Bewegung; die räumliche Bewegung eines 
Massenpunktes; Prinzip von d’Alembert; Energie; 
Impuls; die Planetenbewegung; vorgeschriebene Be- 
dingungen, synthetische Methode; vorgeschriebene 
Bedingungen, analytische Methoden; statisches Gleich- 
gewicht bei Zwangsbedingungen; vorgeschriebene 
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Bedingungen, Lagrangesche Gleichungen zweiter Art; 
Umformung der Lagrangeschen Gleichungen zweiter 
Art; Grundbegriffe der Variationsrechnung; die 
Drehung; beschleunigte Koordinatensysteme. Der 
zweite Teil enthält: Einteilung der Punktsysteme; 
die Kräfte bei einem Punktsystem und das Kraft- 
moment; der Schwerpunkt und der Schwerpunktsatz; 
die Impulssätze; das Impulsmoment und der Flächen- 
satz; das Gleichgewicht der Kräfte; Beispiele zur 
Mechanik der Punktsysteme; der Energiesatz, die 
Lagrangeschen Gleichungen; die Grenzen der klassi- 
schen Mechanik; Grundbegriffe der Festigkeitslehre; 
Gleichgewicht der äußeren Kräfte; die inneren Kräfte 
des Körpers; Kinematik des starren Körpers; Stoß 
von zwei Körpern; das Trägheitsmoment; Grund- 
gesetze der Kreiseltheorie; der kräftefreie Kreisel, die 
Nutation; die Eulerschen Kreiselgleichungen; der 
Kreisel bei äußerem Kraftmoment, die Präzession; 
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die Kreiselgesetze bei drehendem Bezugssystem; An- 
wendung der Kreiselgesetze in Physik und Technik; 
Stoffeigenschaften der Flüssigkeiten und Gase; Statik 
der Flüssigkeiten und Gase; Dynamik der idealen 
Flüssigkeit; Dynamik der zähen Flüssigkeit; der 
Widerstand von Körpern in Flüssigkeiten. de 

Bezüglich des Massenträgheitsmomentes hat der 
Verfasser eine Darstellung für die Komponenten des 
Trägheitstensors gewählt, welche nicht mit einer all- 
gemeineren Tensorrechnung im Einklang steht (bei 
den Zentrifugalmomenten als Tensorkomponenten 
hätte das negative Vorzeichen von vornherein ein- 
geführt werden müssen). - 

Die klare Sprache und die sorgfältige Ausführung 
der Skizzen sind geeignet, dem Leser das Verständnis 
der mechanischen Grundgesetze und ihrer Anwendung 
zu erleichtern. . 


Dresden. H. Neu ber. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


Dr.-Ing. habil. H. Richter, Rohrhydraulik. Ein 
Handbuch zur praktischen Strömungsberechnung. 
Zweite verb. Aufl. XI +328 S. m. 217 Abb., 68 Zahlen- 
tafeln u. 32 praktischen Berechnungsaufgaben. Berlin/ 
Göttingen/Heidelberg 1954. Springer-Verlag. Preis 
geb. 34,50 DM. 


Dr.-Ing. Günther Nahrgang, Zur Theorie des 
vollkommenen und unvollkommenen Brun- 
nens. Aus dem Institut für Hydromechanik, Stau- 
anlagen und Wasserversorgung der T. H. Karlsruhe. 
V-+438. m. 23 Abb. Berlin/Göttingen/Heidelberg 
1954. Springer-Verlag. Preis geh. 7,50 DM. 


Prof. Dr. H. Görtler, Ausbildung und Stellen- 
vermittlung der Diplom-Mathematiker. 
(Herausgegeben im Auftrage der GAMM) 35 S. Mos- 
bach-Baden 1954. Physik-Verlag. 


Carl Forssell, Schubfestigkeit und Schubbe- 
wehrung der: Betonbalken. (Kungl. Tekniska 
Högskolans Handlingar. Transactions of the Royal 
Institute of Technology Stockholm, Sweden, Nr. 78.) 
65. m. 29 Fig., 10 Fotos u. 7 Tafeln. Stockholm 1954. 
Hendrik Lindsträhls Bokhandel i Distribution. Preis 
geh. 10. Kr. 


NACHRICHTEN 


München: Die Kommission für Rechenanlagen der 
Deutschen Forschungsgemeinschaft veranstaltete vom 
26. bis 28. April 1954ein Rundgespräch, dessen wissen- 
schaftliche Vorbereitung, Durchführung und Leitung 
in den Händen von Professor Piloty lag. Das 
Rundgespräch hatte den Zweck, einen möglichst 
engen Kontakt zwischen den Mitarbeitern der von der 
Deutschen Forschungsgemeinschaft geförderten In- 
stitute herzustellen. Insbesondere wurde über aktuelle 
Entwicklungsfragen diskutiert. Die beteiligten Insti- 
tute waren: 

Institutfür Praktische Mathematik der Technischen 
Hochschule Darmstadt (Prof. Dr. A. Walther), 
Max-Planck-Institut für Physik, Göttingen (Prof. Dr. 
L. Biermann), Institut für elektr. Nachrichten- 
technik und Meßtechnik der Technischen Hochschule 
München (Prof. Dr. H. Pilot y), Mathematisches 
Institut der Technischen Hochschule München (Prof. 
Dr.R. Sauer), Institut für Techn. Physik, Universi- 
tät Tübingen (Prof. Dr. Tischner), Institut für 
Theorie der Elektrotechnik der Technischen Hoch- 
schule Stuttgart (Prof. Dr. Bader), und Institut für 
Fernmeldetechn. Geräte und Anlagen der Technischen 
Hochschule Darmstadt (Prof. Dr. Gundlach). 


Zum Rektor der Technischen Hochsehule Mün- 
chen wurde für das Studienjahr 1954/55 der o. Pro- 
fessor für Höhere Mathematik und Analytische Me- 
chanik, Dr. d. techn. Wissenschaften Robert Sauer 
gewählt. 


Dem Privatdozenten für Getriebelehre und Kine- 
matik an der T. H. München, Dr. phil. habil. Rudolf 
Beyer wurde die Amtsbezeichnung ‚„Apl. Professor“ 
verliehen. 


Der wissenschaftliche Assistent am Mathematischen 
Institut der T. H. München, Dr. rer. nat. F. L. Bauer 
wurde zum Privatdozenten für „Reine und Angewand- 
te Mathematik‘‘ ernannt. 


Halle/Saale: In der Nacht vom 8. zum 9. Okto- 
ber verstarb plötzlich an einem Herzschlag der o. Prof. 
der Mathematik an der Universität Halle, Dr. Hein- 
richBrandt, Vizepräsident der Deutschen Akademie 
der Naturforscher, Leopoldina. 


Düsseld’orf: Nach kurzer schwerer Krankheit 
starb am 11.9. 54 Prof. Dr.-Ing. Ulrich Graf. 
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